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ABSTRACT 


A computer program (written in FORTRAN) is presented 
which uses Monte Carlo techniques to simulate one-searcher, 
one-target passive acoustic ASW search that terminates at 
detection. A threshold crossing detection model is used, 
and stochastic variations inthe acoustic Signal are modeled 
uSing either a Lambda-Sigma Jump or Gauss-Markov error 
process. Both platforms have the Sapoeutcrey of detecting 
each other, and area and barrier searches are modeled. 
Features of the program include interactive data input, 
extensive use of graphical displays, and thorough statis- 


tical analysis of the results of the simulation. 
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lee “DESCRIPTEOW OF THE SIMULATION MODEL 


A. INTRODUCTION 


Although the body of analysis devoted to search theory 
is extensive (see, for instance, the lists of references in 
Koopman [Ref. 1] and Washburn [Ref. 2]), most of the tech- 


niques are applicable only to stationary targets. Recently 


there have been some important advances in the theory of 


search for moving targets (see Eagle [Ref. 3], Brown [Ref. 


Stewart [Ref. 5], Stone [Ref. 6], and Washburn [Ref. 7]). 
Even so, it is still a very difficult, if not impossible, 
computational problem to calculate the probability of 
detection or the mean time to detection when the searcher 
has a speed close to that of the target and follows a 


realistic (i.e. not random) track. As a result, computer 


Simulation of the search process is often used to evaluate 


the effectiveness of search tactics. 


4], 


The Monte Carlo simulation computer program presented in 


this thesis is an attempt to model the passive acoustic 


search process in a one-searcher, one-target scenario where 


both the searcher and target platforms have sensors capable 


of detecting each other. The results of the simulation 


allows the analyst to closely approximate parameters (e.g. 


cumulative probability of detection, detection times, detec- 


tion ranges) which are essential in planning search 


operations. With these results, the military planner can 


14 





make informed decisions concerning asset allocation and 
specific platform employment. 

The computer program, Passive Acoustic Search Simulation 
(PASS), 1s based on an acoustic fluctuation model described 
ify MeCabe [Refs 10:pp.7-23] and Hurley [Ref. llspp.25-34]. 
Precedessor programs include the FORTRAN program BEAR, and 
the ALGOL program SEARCH. BEAR was written by L. K. Arnold 
of D. H. Wagner, Associates (see Arnold [Ref. 8] for program 
documentation). SEARCH, a derivitive of BEAR, was written 
by W. J. Browning and J. Risberg, also of Wagner, 

Associates (see Browning [Ref. 9] for program documentation) . 
SEARCH has been used extensively in operational and exercise 
environments under the direction of Submarine Development 
Squadron Twelve (SUBDEVRON-12), where it is currently 
resident on the Semi-Automatic Reconstruction Facility 
(SARF) computer. The basic structure of PASS is similar to 
that of BEAR, with some of the features of SEARCH added. 
Additionally, some modifications were made to algorithms in 
these earlier programs, and new features were added, 
including interactive data input and extensive use of 
graphical presentation of simulation results. 

PASS is written in FORTRAN and is designed for use at 
the Naval Postgraduate School (NPGS). PASS is not portable 
in that it uses external subroutines in the Non-International 
Mathematics and Statistics Library (NONIMSL) and the DISSPLA 


graphics system. The subsequent sections of this chapter 


iS 





are designed to give the reader an overview of the model 
without going into the specifics of the program code or 

logic. If the reader is interested only in learning how to 
use PASS, it 1s recommended that this chapter and Appendix A 
be studied. Appendix A is designed as a User's Guide and 
does not require any knowledge of the program code other 

than variable names, which are described in Appendix B. If 
the reader is interested in the structure of the program and 
the model, then this chapter and Chapter II should be studied. 
Some examples of varied applications of PASS are shown in 


Chapter III. 


B. GENERAL PROGRAM DESCRIPTION 

Program PASS is a Monte Carlo simulation of two passive 
sonar platforms in a search scenario. Arbitrarily, one 
platform is called the "Searcher", and one the "target", 
where the target is typically thought of as a submarine, and 
the searcher aS any paSSive sonar platform. Since both ships 
have passive sensors, the assignment of the names searcher 
and target does not connotate strict roles (1i.e., the searcher 
may be detected first, and becomes a target in the truest 
sense.). 

During each replication, the target’s motion is 
partially random within a user-defined search area. The 
searcher motion is deterministic and controlled by user 


input. Each sensor is subject to random acoustic fluctuations 
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which induce random variations in detection ranges _— 
time. Detection occurs whenever one platform obtains a 
positive signal excess, unless the user employs an integra- 
tion time model, which requires weak Signals to be present 
longer than strong signals for detection to occur. 

The Output of the program consists of distributions of 
parameters of interest (e.g. time of detection and range of 
detection), and point estimates derived from the statistics 
of these distributions (e.g. mean time to detection and 
mean detection range). The format of the output is numerical 


for point estimates, and graphical for distributions. 


fe oeeech SCENARIOS USED IN PASS 
1. Area Search 
The target 1S constrained to move within a specific 
rectangular area ("Search area") in a "semi-random" fashion. 
Target speeds are distributed uniformly over a user speci- 
fied speed range, and courses are uniform on (0, 360) 
degrees. Times between course and speed changes are exponen- 
muamely distributed. The size of the search area and the 
parameters of the random motion are user specified. The 
searcher moves along a track which need not be contained 
within the search area. The Searcher track and speeds are 
user specified. 
2. Barrier Search 
The target is constrained laterally to be within 


certain user specified bounds, thus simulating a "choke point". 


7 








Target motion through the choke point is generally vertical 
(i.e. from “top” to “bottom") with random course and speed 
deviations parameterized by the user. The searcher track 
1S specified as in Area Search. 

The distribution of the initial lateral position of 
the target can be specified by the user as either uniform 
across the choke point, or with Specific probabilities as a 
function of distance across the choke point. In uSing the 
latter option, and by limiting the maximum course deviation, 
the user can control the distribution of the lateral target 


position when the target venetrates the barrier. 


D. THE ENVIRONMENTAL MODEL 

The environment 1S modeled in a conventional manner. 
PASS requires range/propagation loss information for both 
the searcher and target sensors. Direct path (DP) and 
convergence zone (CZ) data are entered Separately. A cubic 
interpolation routine is used to determine propagation loss 
between the user-input data points. The convergence zones 
are modeled as inverted "Square-wells" superimvosed over the 
DP propagation loss curve. Up to five convergence zones are 


allowed. 


fee tHE DETECTION MODEL 
In PASS a detection is always a "Secure detection". That 
is, a detection by the searcher means that the searcher 


detected the target before the target detected the searcher. 
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Simultaneous detections are treated as special cases where 
neither the searcher or the target gets credit for a 
detection. The simulation terminates (replication ends) 
whenever: 

1. The searcher makes a detection. 


2. The target makes a detection (sometimes referred to 
aS a counter-detection). 


3. Awsamultaneous detection occurs. 


4. In the barrier scenario, the target crosses a user 
specified lower boundary. 


5. User specified maximum search time is exceeded. 

The propagation loss curves are used in conjunction with 
searcher and target figure-of-merit (FOM) to determine 
detection ranges. The program uses FOM data which must be 
precalculated by the user based on specific platform source 
levels (Ls), backgroung noise (self noise and ambient noise) 
and sensor directivity (Le), and processor (machine and 
operator) recognition differential (Nrd). 

meocasener Pvqure-Of-Merit 

The Ls of the target 1S assumed to be independent 
of target speed. In effect flow related noise and speed 
dependent propulsion noise are ignored. At the user's 
option, target LS may be made a function of target aspect. 
It is assumed that the deterministic sensor performance of 
the searcher is affected by the searcher speed. This results 
in a searcher FOM that is dependent on Searcher speed, and, 


optionally, target aspect. 
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Ze Target pr gaure=o f=Marit 


The Ls of the searcher 1S assumed to be speed depen- 
dent, and independent of searcher asvect. This results ina 


target FOM that is dependent on searcher speed. 


3. The Threshold Crossing Model With a Stochastic Error 


Process 

Given a FOM against a target for a search platform, 
the maximum detection range of the searcher is determined 
fromthe provagation loss curves. This results in the sensor 
conforming to the "definite range law", or "cookie-cutter" 
performance. This type of sensor readily lends itself to 
mathematical analysis, but has the disadvantage of not 
reflecting actual sensor performance. What 1s seen under 
operational conditions 1S a variation in Signal strength 
(and thus a variation in detection ranges) over time, even 
fHOretargets at fixed range and source level. 

To model these variations in received acoustic 
Signal, a stochastic error process is added to the determin- 
istic figure-of-merit. This results in the signal excess 
consisting of deterministic and random components. For a 
detailed discussion of threshold crossing models and 
associated stochastic error processes, see McCabe [Ref. 10] 
and Hurley [Ref. 11]. See Appendix D for a discussion of 
the relationship of signal excess (SE) to figure-of-merit 
(FOM) . 

The mean signal excess (SE) is defined as the 


expected difference between the signal-to-noise ratio (SNR) 
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at the hydrophone output (in decibels, db) and the SNR that 
is estimated to afford a probability of detection of 0.5. 
The model used for the "observed" signal excess, as 


e “fanction of tame in PASS, is of the form: 


SE = SE + xX+ Y (ea) 


where SE is the deterministic mean Signal excess, and X and 
Y are random variables that describe the fluctuations in the 
Signal level. 

If we let the subscript s apply to the searcher and 


t for target, the model used in PASS is: 


SE_(t) = SE,(t) + X(t) + Y(t) G2) 


SE, (t) = SE, (t) + ¥(t) + Z(t) ies) 


where: 


1. X(t) is a stochastic error process which describes the 
fluctuations "local" to the searcher. 


2. Z(t) is a stochastic error process which describes the 
fluctuations "local" to the target. 


3. Y(t) is a stochastic error process which describes 
acoustic fluctuations of a global nature whose effect 
1S common to the searcher and target. 

4. X(t), Y(t), and Z(t) are independent processes. 

X(t) and Z(t) may be thought of as modeling the 
fluctuations that affect only the performance of the 


detecting platform, such as onboard fluctuations in sensor/ 


processor performance, including operator performance. 


Z1 








Y(t) may be thought of as modeling the gross 
acoustic fluctuations that occur along the “environmental 
path" between the platforms, such as temperature and 
Salinity gradients and changes in bottom conditions. 

With signal excess now a function of the mean signal 
excess and the stochastic effor processes, we have the means 
to model the behavior of the sensor by employing a threshold 
crossing detection model as follows: 


1. If SE(t) is less than 0, the probability of detection 
#5 0.0. 


ey wt SE(t) iS greater than or equal to 0, the probability 
Of detection is 1.0. 


tas, a G@etection occurs at time t if, and only if, SE(t) is 
greater than or equal to zero. The fact that SE(t) is now 
undergoing random fluctuations about the mean Signal excess 


results in random variations in tne sensor detection range. 


meee a LOCHASTIC ERROR PROCESSES 

PASS allows the user to choose either a Lambda-Sigma 
Jump (LSJ) or a Gauss-Markov (GMA) error process to model 
the acoustic fluctuations. Hurley [Ref. 11] gives an over- 
view of these processes, andothers used in search simulation 
programs. McCabe [Ref. 10] provides a more mathematically 
complete treatment of the processes, along with a compre- 
hensive comparison of the LSJ and GMA models. In PASS we 
use the sum of two LSJ or GMA process, frequently referred 
to aS a compound process. In the following discussion, a 
simple process, vice the compound process is considered for 


clarity. 
22 





l. The Lambda-Sigma Jump Error Process 


Define a Lambda-Sigma Jump error process, &(t), as 


follows: 
Eers) = 20%) 6(t)+[(Lez(t) )n Crud) 
(ees T ) 
i me <7 
z(t) = (1945 ) 
Oita = T 
where: 
1. t 1S an exponential random variable with rate parameter \. 
2. ni1S anormal random variable with zero mean and 


VaGiance o- 

The process can be though of as a "marked" Poisson 
process with rate parameter \, where the magnitude of the 
marks are themselves random variables with a normal distri- 
ipmeon (with zero mean and variance o*). The standard 
deviation of the normal distribution, o, can be thought of 
aS a scale parameter in this process. The figure-of-merit, 
then, 1S constant over exponentially spaced time intervals, 
and changes to new levels determined by a normal distribution 
at the end of each of these time intervals. 


Themeenarcmance function for €(t), 1s: 


Couratemn ec (crs) |= o2e™ > (1.5) 


This function iS a measure of the correlation between values 
of &€(t) at different times. From this we can see that §&(t) 


is second order stationary. That is, the covariance depends 
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on the time difference, s, and not the time, t. The para- 
meter A determines the amount of dependence between successive 
acoustic levels, and can be chosen to yield any result from 
independence to complete dependence. The resulting process 
1s then a piecewise continuous function in which there are 
periods of complete dependence of detection opportunity 
interrupted by fluctuations that introduce independence 
between these periods. 

In PASS, then, we must provide a total of six param- 
eters to completely describe the compound error processes 
used to model the acoustic fluctuations. They are: 


Ae eectemDasamercr fomgrhe searcher-—local error 
process, X(t). 


eo Seenespanateter for the searcher—local error 
process, X(t). 


ee eeesace Patameter for Ehe global error process, Y(t) . 
4. g, = scale parameter for the global error process, Y(t). 
5. A,= rate parameter for the target local error 


3 
process, 42(t). 


6. o, =scale parameter for the target local error 
process, Z(t). 


The selection of appropriate values for X} and o isa 
not well understood function of the environment, sensor, and 
processor. In practice, the selection of the values are 
subjective decisions based on the experience of the analysist 
as what seems to yield reasonable results. Tehan [Ref. 12] 


provides a discussion of estimation techniques for 
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these parameters. In applying these recommendations to 
PASS, which has compound processes, the values of A and o 


were governed by: 


2 
BSA dyes 6S(a{+02) <9 
(lieaGy) ~ (1.7) 
3S +r <5 6S(0o +o ) S9 
Z 3 2 3 


See Appendix E for a description of the computer 
Simulation of LSJ. McCabe [Ref. 10:pp.9-13] gives some 
analytical results based on the LSJ process. 

2. The Gauss-Markov Error Process 

The Gauss-Markov process, u(t), introduces explicit 
dependence of the value of u(t+s) based on the value of u(t) 
for any s. The process is Markovian in that the present 
value depends only on the value that immediately proceeds 
it. The process is GausSian in that, for any n, the joint 
eee outiOn Of {&(ti)}, i=l,...-,n is multivariate-normal 

2 


each with zero mean and variance o . In Simulating the 


process, we can use the functional form: 


u(tts) = @ ASu(t) + kn (1.8) 
where: l 
Z ay 
meet wer (see Appendix E) 
2. mn is a normal random variable with zero mean and 


varlance o 
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The process results in a continuous sample path, but 
6, in fact, nowhere differentiable. The desirability of 
using this type of error function in modeling acoustic fluc- 
tuations is that it 1S "more like reality". Whether it is in 
any sense "better" than the LSJ process has not been 
determined. Results of other simulations indicate that 
cumulative probabilities of detection are higher, and mean 
time to detect is lower, when the GMA process 1S used instead 
of the LSJ. However, the differences are not Significant to 
the point of determining which process 1S most appropriate 
aS an acoustic model. One drawback of the model is that it 
1s not possible to exactly duplicate a continuous function 
on a digital computer, and thus the machine representation 
of the GMA is only a close approximation of the continuous 
error process. Thus, the manner in which the GMA process is 
digitally replicated has an effect on the results of the 
Simulation. 


The parameterization of the process is exactly 


analogous to that of the LSJ process. Again, N yg can be 
a f 
thought of as rate parameters, and 9, , can be thought of 
f f 
aS scale parameters, exactly as in the LSJ model. In some 


analyses, 1/A is referred to as the “relaxation time" (and 
A as the "relaxation coefficient") since it measures the 
decay time of statistical dependence between random variables 


in the fluctuation process for both the GMA and LSJ processes. 
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The covariance function for the Gauss-Markov process 
is exactly the same as that for the LSJ (see Equation 1.5), 
and therefore, the GMA is also a second order stationary 
process. 

There are no analytical results similar to those 
from the LSJ process when the Gauss-Markov process is used. 
See Appendix E for a description of the computer representa- 


tion of the GMA process. 


feo a Peerets—-OUT—OF-FIVE DETECTION CRITERIA MODEL 

McCabe [Ref. 10:p.50] contends that when a target is 
exposed for only a short amount of time but at a relatively 
high signal-to-noise ratio (SNR), the threshold crossing model 
returns an unrealistically high probability of detection, and 
presents a signal integration model based on a time dependent 
recognition differential. PASS has, as an option, a much | 
simplified model which is easier, and faster to implement. 

It seems necessary to alter the threshold crossing model 
under the above conditions so that detection would require 
some minimum combination of signal strength and time. This, 
in effect, would model what is actually integration time for 
the sonar processor. In general, the desired relationship 
is that the stronger a Signal is, the shorter the time 
required for recognition as a valid contact. This model 
seems to be especially necessary when detection is based on 


other than aural cues, as is the case in modern sonars. 
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As a user option, instead of the threshold crossing 
model (where detection occurs the first time SE(t) is equal 
to or greater than zero), a 3-of-5 model can be employed 
wich sequires SE(t) = 0 on 3 of the last 5 samples of 
Signal level for detection. This Minimum Signal Excess 
Logic Model (MSEL) results in the desired effect of a weak 
Signal being required to be "present" for a longer period 
for detection to occur. See Appendix F for a more complete 
analysis of the model. 


McCabe [Ref. 10:p.50] recommends that an integration 


model be used whenever strong convergence zones are present. 


MH. "A SEARCH MEASURE OF EFFECTIVENESS (MOE) 
1. Searcher MOE 
We may wish to measure the effectiveness of search 
(Or evasion) in terms other than the probability of secure 
detection. For example, we may wish to penalize the searcher 
heavily for a secure target counter-detection, less heavily 
for a simultaneous detection, and still less heavily for no 
detection (escape). 
nee 
1. Ns = the number of secure searcher detections 
2. Nt = number of secure target detections 
3. Nb = the number of simultaneous detections 
4. Nn = number of no detection replications 


5. Nr = the total number of simulation replications 
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Then define a MOE for the Searcher as: 


_ Ns~3Nt-2Nb-Nm 
MOES = —S CC i Cee > (fi 39) 


where 
OSMOESS1 and MOEs = 1 => Ns = Nr 


Note that the weighting factors are a function of the 
relative penalties assigned to Situations other than a 
secure detection, and could be easily adjusted. 
2. Target MOE 
In the area search scenario, we assume the target to 
be aggresSive, and therefore a target MOE 1S assigned analo- 
gous to the searcher MOE: 
Nt-3Ns-2Nb-Nm 


In the barrier search, the target 1S assumed to want 
to avoid contact with the searcher, so the target MOE is 


adjusted as follows: 


_ Nt-3Ns-2Nb+Nm 
MOE | ——- 1.45 teOh5 (ee) 


where, in both cases, 


OSMOE, Sl and MOE, =~ l=?N, =Nr 
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3. The Exchange Ratio 


The exchange ratio is defined as: 


MOE. 


Which provides a composite MOE for the search. 


In area search: 


ER>Ll = NS oNt 

ER = 1 Ns=Nt 
ER 1 = Ns<Nt 

In barrier search: 


il 


ER>l => Ns >Nt+5Nn 
BR=- l= Ns=Nt+5Nn 


ER<l <=> Ns<Nt+5Nn 


I. PROGRAM RESULTS/OUTPUT AND STATISTICAL ANALYSIS 
The following resultS are available from PASS: 
1. A complete record of the input data. 
2. Probability of detection by the searcher 
a. Fraction of detections wnich were direct path. 


b. Fraction of detections which were in the 
convergence zones. 


c. Fraction of detections in each convergence zone. 
3. Probability of counter-detection by the target. 
a. Fraction of counter-detections which were DP. 


b. Fraction of counter-detections which were in the 
CZs. 
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e- Bea@etion of courteradeteCtioms in each GZ. 
Histogram and statistics of time to detection by searcher. 


Histogram and statistics of time to counter-detection 
by the target. 


Histogram and statistics of range at detection by 
searcher. 


Histogram and statistics of range at counter-detection 
by target. 


Sectioning of time and range data (if sample size is 
large enough) to allow for estimates of the variability 
of sample statistics. 


A plot of cumulative probability of detection versus time. 


A plot of cumulative probability of counter-detection 
versus time. 


A plot of conditional probability of detection versus 
range. 


A plot of the searcher positions when the searcher 
detects the target. 


A plot of the target positions when the searcher 
detects the target. 


Perot Of CONG tional probability of “counter=—deteetion 
versus range. 


A plot of target positions when target counter-detects 
searcher. 


A plot of the searcher positions when the target 
counter-detects the searcher. 


A plot of typical searcher and target FOM for a five 
hour period. 


J. VALIDATION OF THE PROGRAM 


There is no existing body of data or analysis that would 


lend itself to a direct comparison to PASS results for a 


rigorous validity check. Numerous area and barrier search 
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Scenarios, with various environmental, platform, and 
"random" parameters, were run with results not inconsistent 
with operational experience, limiting theoretical results, 
and other simulation program results. Due to the large 
number of variations the program allows, not all could be 
tested. The strongest statement that can be made regarding 
PASS validity is that it has not shown to be invalid after 


running in various configurations. 


K. PROGRAM RUN TIME 

PASS is a relatively expensive analysis tool from the 
Standpoint of computer run time. On the NPGS computer (IBM 
VM-370 with three 3033 CPU's) a rule of thumb is that it 
will take about 50 CPU-seconds per 1000 replications of an 
area search in 100 nm x 100 nm Search ec. Obviously, more 
or less time will be required for different sensors and 
other scenarios. The point to be made is that the run time 
is long enough to warrant extenSive pre-planning so that the 


Meeimum desired information is obtained from each run. 
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II. DESCRIPTION OF PROGRAM PASS 


Ae INTRODUCTION 

The basic flow diagram for program PASS is shown in 
Poammees 2.1, 2.2, and 2.3. Each “part" of the program, as 
described in the comment blocks in the source code (see 
Appendix H) iS indicated in parentheses in each block of the 
flow diagram. The flow diagram, the annotated source code, 
and the description of the variable names (See Appendix B) 
Should be used with this chapter to gain an understanding of 
the logic employed in the simulation and the imolementation 
of this logic in the FORTRAN code. 

B. DATA INPUT AND SIMULATION INITIALIZATION (PROGRAM 

Permts 1, 2, 3) 

All logical variables are set to default values prior to 
data input. 

Data input is normally accomplished by the user respond- 
ing to a number of ovtions presented on the terminal screen. 
These options are displayed by sequential calls to sub- 
routines OPTNO through OPTN20. During this interactive 
session, the data necessary to run the program in the 
default mode is written to file 07 on the A-disk, which 
allows saving the basic data from one run to the next. The 
user can bypass all the ovtions on the screen and proceed 


directly to program execution by reading the data from this 
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Bebe through the subroutine READIT. Invoking subroutine 
READIT requires accepting the following defaults: 
= = complete, did Dbroperly formatted, inpwt data set is 
assumed to be on the A-disk in file PASS DATA (file 
Gy). . 


2. All input data will be sent to the output file PASS 
C@PrUTY (tite 6) . 


3. The acoustic fluctuation model is Lambda-Sigma Jump 


iss). 

4. A complete statistical analysis of the results is to 
be done. 

>. the 3=-of=-5 detection criteria model is not used. 

6. For area search the initial target position is 


distributed uniform on (0,XMAX), (0,YMAX). 


7. For barrier search, the initial target lateral position 
1s distributed uniform on (0, XMAX). 


8. The searcher figure-of-merit 1S independent of target 
eaaeecct. 


Subroutines ECHO] through ECHO6 write the input data to 
the output file if the user so desires. 

The limits of the target course variations in the 
barrier scenario (ANG1, ANG2) are initialized if barrier 
search is selected. Detection counters (NDO, NDT, NCZDO(I), 
NCZDT(I), NBOTH, NONE) are set to zero, and the stack pointers 
(Cl, C2, C3, C4) are initialized. The direction of each 
Searcher path leg (DX(I), DY(TI)) and the distance of each 
-leg (DIST(I)) are computed, anc the target speed range 
(STINC) and boundary reflection constants (TOXMAX, TOYMAX) 


are calculated. 
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C. REPLICATION COUNTER AND SETUP (PROGRAM PART 4) 

If the required number of replications (NREP) have been 
completed, transfer is made to the data analysis and output 
routines (Part 18) and subsequent program termination. 
Otherwise, another replication is started. 

A count of completed replications is written to the 
screen every 200 replications. Counters for the collection 
of a representative sample of figure-of-merit data (Ml, M2) 
are reset if five hours of this data has not been collected 
in a previous replication. Simulation times (TNOW, TLAST) 
are set to zero, and signal MSEL counters (MSELO, MSELT) and 
convergence zone (CZ) detection pointers (KC2Z0, KCZ2T) are 


reset. 


D. INITIALIZATION OF THE REPLICATION (PROGRAM PART 5) 
At the start of each replication, the following initial- 
ization and setup is accomplished: 


1. The time to the next fluctuation change (TIFL(I), 
T=1,2,3) is selected from an exponential distribution, 
parameterized by ALAM(I), by calls to subroutine EXPO. 


2. The initial magnitude of the fluctuation (AFL(I)) is 
selected from a normal distribution, parameterized by 
SIGMA(I), by calls to the subroutine XLS. 


3. The searcher position (XO, YO) is set to the first 
searcher track anchor point (XP(l1), YP(1)). 


4, The search leg pointer (NLEG) 1s set to one, and the 
search mode (MODE) is set to "drift". The time to the 
next searcher speed change (TSC) is set to the drift 
time. The figure-of-merit for both platforms (FOMO, 
FOMT) is) selected based on searcher drift speed, and 
searcher speed (SO) is set to drift speed. 
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The time to the searcher course change (TCC) is set 
based on the first leg distance and drift speed, and 
Searcher speed vectors (VX, VY) are set based on the 
direction vectors and search speed. 


If the user is to specify the initial target lateral 
position in the barrier scenario, the initial target 

X-position (XT) is selected by a call to subroutine 

XDISTB. Otherwise, the initial target X-position is 

selected uniform on (0,XMAX) by a call to subroutine 
UZ. . 


If the barrier scenario is selected, the initial 
target Y-position (YT) is selected a distance offset 
from YMAX. This offset distance is representative of 
the searcher being "“time~-late" on the barrier. The 
offset is calculated based on the maximum time late 
(START), the mean target speed, and a call to sub- 
routine UZl. In the area search scenario, the initial 
target Y-posSition is selected uniform on (0,YMAX) by 

a call to subroutine UZ1l. 


Target speed (ST) is selected uniform on (STMIN,STMAX) 
by a call to subroutine UZ1l. 


The time to target speed change (TTSC) is selected 
from the exponential distribution, parameterized by 
RESC, by a call to subroutine EXPO. 


In the barrier scenario, the initial target psuedo- 
course (X) is selected uniform on (ANG1,ANG2), and in 
the area search scenario, uniform on (0,360), by a 
call to subroutine UZl. If the searcher FOM is 
dependent on target aspect, the target psuedo-course is 
saved (THETA). Target direction vectors (COSX, SINX) 
and speed vectors (UX, UY) are computed based on 
target course and speed. 


The time to target course change (TTCC) is selected 
from an exponential distribution, parameterized by 
REee, bya Gall to subroutine EXPO. 


DETERMINE DETECTION RANGES (PROGRAM PART 6) 


If the searcher FOM is dependent on target aspect, the 


searcher FOM computed in Part 5 or 14 is re-calculated based 


on the current relative bearing from the target to the 


searcher (BREL). The subroutine RELB returs BREL based on 
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the current searcher-target geometry. The subroutine INTRPL 


returns interpolated values of searcher FOM (FOMO) based on 
BREL and the searcher-speed/target-relative-bearing FOM data 
(FOMBD(I), FOMBS(I)). The searcher FOM (X) is computed by 
applying the error process values (AFL(1), AFL(2)) to FOMO. 
Five hours of searcher FOM (FO(I)) is saved for graphical 
display. The direct-path (DP) detection range for the 
searcher (RNGO) is obtained by applying the searcher FOM (x) 
to the DP propagation loss curve through subroutine INTRPL. 
If convergence zones are present (NCZO not zero), the 
searcher FOM (X) is applied to each CZ propagation loss value 
(CZLO(I)) to determine the most distant CZ in which a detec- 
tion can take place (KCZO). If KCZO is zero, then no 
Searcher detection can take place ina CZ. 

The preceeding procedure, with the exception of the 
aspect dependency portion, is applied also to the target 
related data to produce values for RNGT and KCZ2T. 

The maximum range at which a direct-path detection can 
take place by either platform (RMAX) 1s calculated from the 


maximum of RNGO, RNGT. 


F. CHECK FOR DETECTION (PROGRAM PART 7) 

The range between the platforms (RNG) is calculated 
based on the current geometry, and the algebraic difference 
between this range and the maximum DP detection range is 
Calculated (TBIG). Direct-path detection is possible for at 
least one platform only if TBIG is less than or equal to Zero. 
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If the searcher FOM will support detection in one or 
more CZs (KCZO not zero), the range (RNG) is checked to see 
it the target is in amy searcher CZ annulus, starting from 
the most distant one that will support detection, and 
working in toward the searcher. If the target is in one of 
the CZs, KCZO is set to minus-one (as a CZ detection flag), 
and the CZ in which detection 1s possible (ICZO) is saved. 

An identical CZ detection check for the target is done, 
and if the target can detect the searcher in a CZ, KCZT is 
set to minus one, and ICZT is saved. 

The value of T8SIG is set to the minimum distance the plat- 
forms would have to traverse (asSuming a head-to-head closure) 
for a detection to just occur, based on the current searcher 
and target speeds and FOMs. That is, if the searcher and 
target were headed directly toward each other, TBIG would be 
the total distance covered by both platforms at the time of 
meeete detection by either. 

If a CZ detection iS possible (KCZO and/or KCZT less 
than zero) or if a DP detection is possible (TBIG less than 
zero), a transfer is made to Part 16A where it is determined 
1f a detection 1S made, and by whom. 

If detection is not possible, transfer 1S made to 


Part 8. 


G. COMPUTE THE TIME OF THE "NEXT EVENT" (PROGRAM PART 8) 
If the MSEL model is used, and either MSEL counter (MSELO, 


MSELT) is not in a zero state, this implies that detection 
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was possible at some earlier time, and i int 2@ enw. 
Therefore, any non-zero MSEL counter is decremented by one. 
If, after decrementing, any MSEL counter is not in a zero 
State, TBIG is set to zero. This results in the maximum 
Simulation time-step defaulting to 0.05 hours. 

The maximum simulation time-step is calculated by taking 
the maximum of 0.05 hours (3 minutes) or the time required 
for searcher and target to close the distance TBIG assuming 
all of the current platform speeds (SO, ST) are in the line 
of sight (head-to-head closure). The range TBIG is replaced 
by the time TBIG. The incremented time (TINC) 1S set to the 
current time (TNOW) plus TBIG, and the next event pointer 
(J) 1S set to 8. Next, the times of all other possible 
events are checked to see if any are earlier than TINC. If 
an event time 1s earlier than TINC, the time of this event 
replaces the current value of TINC, and the next event 
pointer is changed appropriately. When this process is 
completed, TINC will be the time of the "next event", and 
the value of J will identify what that event is. Note that 
if J remains at 8, the time increment will be TBIG. TBIG 
tends to be large whenthe platforms are far apart, and 
decreases as the range decreases. 

If TINC exceeds the maximum allowable search time (TMAX), 
then transfer is made to Part 17 for the recording of no- 
detection data, and the start of a new replication, if 


appropriate. 
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H. MOVE THE SHIPS (PROGRAM PARTS 9A, 9B) 

The time increment (X) is obtained by subtracting the 
current time from the incremented time, and the projected 
searcher position (XOT, YOT) is obtained by applying the 
time increment to the searcher velocity vectors (VX, VY). 
The distance from the current searcher position (XO, YO) to 
the projected searcher position (Dl) is calculated from the 
geometry. The distance fromthe current searcher position 
to the next searcher track anchor point (D2) is calculated 
from the geometry. If Dl is less than or equal to D2, the 
searcher position is updated to the projected posmtwen. iIf 
Dl is greater than D2, then the time increment (X) is adjusted 
to put the searcher at the next anchor point, and the next 
event pointer (J) is set to 4 which will, in Part 12, cause 
the searcher position to be updated to the anchor point. 

The target position (XT, YT) 1S updated by applying the 
time increment (X) to the target velocity vectors (UX, UY). 

The simulation time (TNOW) is updated by the time 
increment (xX). 

If the target position, as calculated in Part 9A, ends 
up outside the search area (0,XMAX), (0,YMAX), the target 
Peerelon and velocity vectors are adjusted such that the 
boundaries appear to "reflect" the target. The target track 
behaves as would a light beam striking the mirrored surfaces 


of the search area boundaries. If the searcher FOM is target 
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aspect dependent, then the new target course must be saved 
(THETA) for use in subsequent calculation of relative 


bearing. 


I. BRANCH TO THE "NEXT EVENT" (PROGRAM PART 10) 

Depending on the valueof the next event pointer (J), 
transfer 1S made to the program part which accomplishes the 
physical event associated with the time-step computed in 
Pere 3. 

If J=8 (which means "no event" takes place), transfer is 
made to Part 6 or Part 7, depending on the integration model 
in effect. If the MSEL (3-of-5) model is in effect, 
transfer must be made to Part 6 to reset the KCZO and KCZT 
flags. In future sections where transfer 1S made to Part 6 
or 7, the same dependence on the MSEL detection model in 
effect will govern where transfer 1S made. 

i J=1,2,3 then an acoustic fluctuation level is to 
change, and transfer is made to Part ll. 

If J=4 then the searcher course is to change at a 
Searcher track anchor point, and transfer 1S made to Part l2. 

If J=5 then the target course is to change, and transfer 
1s made to Part 13. 

If J=6 then the target speed is to change, and transfer 
toemeace to Part 15. 

Tf J=7 then the searcher speed is to change, and 


Pmemeter 1S made to Part 14. 
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J. ACOUSTIC FLUCTUATION LEVEL CHANGE (PROGRAM PART 11) 

The fluctuation level change is calculated differently 
depending on the model used. The next event pointer (J) also 
identifies which error term is to change. 

Jo tne LoJ medel, a call to the subroutane XLS with 
input parameter SIGMA(J) returns a new fluctuation level 
(AFL(J)) from a normal distribution. Next, a call to the 
Subroutine EXPO with input parameter ALAM(J) returns a new 
time increment (T) to the next fluctuation, which is added 
to current time to obtain the time ofthe next fluctuation 
Oeertis)). Transfer is then made to Part 6. 

In the GMA model, the time to ane next fluctuation 1s 
obtained in a manner identical to the LSJ model. The lag 
time from the last calculation of error signals (S) 1s 
obtained by subtracting the last time the levels were calcu- 
lated (TLAST) from the current time, and then setting TLAST 
to the current time. All three of the signal fluctuations 
are calculated by successive calls to subroutine XLS parame- 
terized by SIGMA(JJ), JJ=1,2,3, and the algebraic expression 
developed in Appendix E which involves ALAM(JJ), JJ=1,2,3. 


Transfer is then made to Part 6. 


K. SEARCHER COURSE CHANGE (PROGRAM PART 12) 

The searcher position (XO, YO) is assigned to the next 
search leg start point (XP(NLEG), YP(NLEG)) and new searcher 
velocity vectors (VX, VY) are calculated based on the next 


leg direction vectors (DX(NLEG), DY(NLEG)) and current 
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search speed (SO). The time to the next Searcher course 
change (TCC) is calculated based on the leg distance 
(DIST(NLEG)) and current searcher speed. Transfer is made 


to Part 6 or 7, depending on the detection model in effect. 


L. TARGET COURSE CHANGE (PROGRAM PART 13) 

If the area search scenario 1s in effect, the new target 
course (X) is selected from a uniform distribution on 
(0,360) by a call to subroutine UZ1l. If the barrier scenario 
is in effect, the new target course 1S selected froma 
uniform distribution on (ANG1,ANG2). The new target seueae 
1s saved (THETA) if the searcher FOM is dependent on target 
aspect. 

New target speed vectors (UX, UY) are calculated based 
on the new course and current target speed (ST). The time 
interval tothe next target course change (T) is obtained 
from an exponential distribution by calling the subroutine 
EXPO with input parameter RTCC, and the time of the next 
course change (TTCC) is calculated based on this interval. 

If the searcher FOM is a function of target aspect, 
transfer is made to Part 6 (via Part ll to calculate new 
fluctuation levels if GMA is the fluctuation model) to 
determine new detection ranges. Otherwise, transfer 1s 
made to Part 6 or 7 depending on the detection model in 


effect. 
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M. SEARCHER SPEED CHANGE (PROGRAM PART 14) 

The MODE state is changed to reflect the new searcher 
$peed. The Searcher and target FOM (FOMO, FOMT) are changed 
to reflect the new searcher speed, and the search speed (SO) 
1s updated. The time of the next searcher speed change (TSC) 
1s updated using the equivalence variable TIME(MODE). The 
time to the next searcher course change is updated based on 
the new search speed, and distance to the next searcher 
traek anchor point (D). Transfer is made to Part 6 (via 
Part 11 if GMA is the fluctuation model) to determine new 


detection ranges. 


N. TARGET SPEED CHANGE (PROGRAM PART 15) 

The new target speed (ST) is selected from a uniform 
distribution on (STMIN, STMAX) by a call to the subroutine 
UZ1. The time interval to the next speed change (T) is 
selected from an exponential distribution by a call to the 
subroutine EXPO with input parameter RTSC, and the time of 
the next speed change (TTSC) is calculated based on this 
interval. New target speed vectors (UX, UY) are calculated 
based on the new speed and current target direction vectors 
(@@sx%, SINX). Transfer is then made to Part 6 or 7, 


depending on the detection model in effect. 


O. DETERMINE WHICH PLATFORM CAN DETECT (PROGRAM PART 16A) 
This part of the program is entered only if, in Part 7, 


detection by at least one of the platforms is possible. 
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Therefore, if the direct path detection range of the target 
(RNGT) is less than the current range (RNG), and the target 
cannot make a CZ detection on the searcher (KCZT greater or 
equal to zero), the detection opportunity must exclusively 
belong to the searcher, and transfer is made to Part 16C. 
Similarly, if the direct path detection range of the 
searcher (RNGO) is less than the current range, and the 
searcher cannot make a CZ detection on the target (KCZO 
greater or equal to zero), the detection opportunity must 
exclusively belong to the target, and transfer is made to 
Past 16D. 

If neither of the above cases are true, then the detec- 


tion must be simultaneous, and transfer is made to Part 16B. 


P. SIMULTANEOUS DETECTION OPPORTUNITY (PROGRAM PART 16B) 
If MSEL is not in effect, the simultaneous detection 
counter (NBOTH) is incremented, and transfer is made to 
Part 4 to start a new replication or output results. 
If the MSEL model is used, both the searcher and 
target counters (MSELO, MSELT) are incremented. If both 
the counters equal three, then the actions in the preceeding 
paragraph are carried out. If only the searcher counter is 
at three, then transfer is made to Part 16C. If only the 
target counter is at three, transfer is made to Part 16D. 
If no. detection takes place (both counters less than 3) 
transfer is made to Part 8 to determine the time-step to 


the next event. 
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SEARCHER SECURE DETECTION OPPORTUNITY (PROGRAM PART 16C) 


If the MSEL model is not used, the following takes place: 


Searcher detection counter (NDO) is incremented. 

If this was a CZ detection (KCZO = -1), the number of 
detections in the particular CZ is incremented 

GIGZDO (ICZO)). 

The time of the detection is saved (NTDO(NDO)). 

The range of the detection is saved (RNTDO(NDO)). 


Searcher position at time of detection is saved 
(XODT(NDO), YODT(NDO)). 


Target position at time of detection is saved 
iFrODI(NDO), YTORDM( NBO) ) . 


Transfer is then made to Part 4 to begin a new 
Peplication of output results. 


If the MSEL model is used, the following takes place: 


Searcher MSEL counter (MSELO) is incremented. 


If MSELO is 3, the actions in the preceeding paragraph 
are executed. 


If MSELO is less than three, the target MSEL counter 
(MSELT) is decremented if it was not zero, and 
transfer 1S made to Part 8. 


TARGET SECURE DETECTION OPPORTUNITY (PROGRAM PART 16D) 


If the MSEL model is not used, the following takes place: 


Target detection counter (NDT) 1S incremented. 

If this was a CZ detection (KCZT = -l), the number of 
detections in the particular CZ is incremented 

(MCTET (WEZT)). 

The time of the detection is saved (NTDT(NDT)). 

The range of the detection is saved (RNTDT(NDT)). 


Searcher position at time of detection 1s saved 
ODO (NDT) ,@=YTDO(NDT)) . 
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5. 





Target position at tame of detection is saved 
(RETO (NST) , YOTD (NBT) ). 


Transfer is then made to Part 4 to begin a new 
replication or output results. 


If the MSEL model is used, the following takes place: 


Target MSEL counter (MSELT) is incremented. 


If MSELT is 3, the actions in the preceeding para- 
graph are executed. 


If MSELT is less than three, the searcher MSEL 
counter (MSELO) is decremented if it was not zero, 
and transfer 1s made to Part 8. 


NO DETECTION OCCURS (PROGRAM PART 17) 


If no detection occurs by the maximum allowed search 


time (TMAX), the no-detection counter (NONE) is incremented, 


and transfer is made to Part 4 to start a new replication 


Or output results. 


ee 


RESULTS OUTPUT (PROGRAM PART 18) 


By calling subroutine SINKEM, the following output 


functions are accomplished: 


dl 


Some results of the simulation (times, ranges) are 
written to file 06 on the A-disk in the form of 
histograms and statistics. 


Some results of the simulation (times, ranges, 
cumulative probability of detection/counter- 
detection, geographical positions, and figures-of- 
merit) are written to file 10 on the A-disk. File 10 
is designed to be a data file for the FORTRAN program 
PASPLT which produces graphical output by invoking 
the DISSPLA graphics system. 


An abbreviated summary of the Simulation results 1s 


sent to the terminal screen (file 08) and to file 06 
on the A-disk. 
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Figure 2.1. Flow Diagram for Program PASS Parts 1-7. 
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Ill. EXAMPLE SEARCH PROBLEMS 


=. INPRODUCTION 

In this chapter a set of search scenarios will be inves- 
tigated to demonstrate the operation and capabilities of 
PASS. In each scenario (with the exception of the definite 
range law examples), environmental and performance parameters 
are representative of "“real-worid" conditions, but they do 
not represent any specific set of platforms or environmental 
conditions. This was done to preclude questions regarding 
security classification. 

The vurpose of this chapter 1s demonstration. As such, 
the problems investigated are treated in what may be called 
a superficial manner. That 1s, the questions posed are not 
rigorously analyzed, and conclusions drawn can only be 
applied within the narrow constructs of the specific scenario. 

With one exception, the output from PASS (and PASPLT) as 
Shown in the figures in this chapter 1S unaltered. The 
exception is the histogram and statistics output produced by 
the NONIMSL subroutine HISTGP, which 1s produced on the 
computer center line printer. The format of the histograms 
and statistical output 1S oversize, and is not of high 
enough quality to be photo-reduced for legible inclusion in 
this thesis. Selected statistical results from HISTGP will 


be provided in tabular form when required. 
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B. AREA SEARCH: NOISE QUIETING VS. IMPROVED SENSOR 
PERFORMANCE 


1. Se@mario Description 

A target with a relatively large acoustic disadvan- 
tage is the subject of systematic area search. The question 
examined is which of the following possibilities would be 
more advantageous to the target: 

1. Decreasing target radiated noise by 3 db, or 
2. Increasing the FOM of the target Sonar by 3 db. 
= @2ASS Input Data 

Table I shows the input data for the base case (case 
A-I), in the format presented in the PASS output file. 

3. Environment and Search Plan 

Figures 3.1 and 3.2 show the propagation loss models 
used in the example. Note that even before figure-of-merit 
is calculated, the searcher has an advantageous Situation in 
that his propagation loss curve predicts longer detection 
ranges than for the target (for the same FOM). This occurs 
because the geometric mean frequency of the searcher sonar 
is less than that of the target. Due to the presense of 
convergence zones, the MSEL model is used. 

Figure 3.3 shows the searcher track, with track 
spacing of 75 miles (about equal to the range of the second 
CZ). The search area is 300 nm square. 

4. Platform Parameters 
The search speed is established as constant at 10 


knots by setting sprint and drift speeds equal to 10 knots, 


m2 





and searcher sprint and drift FOM equal at 87 db. Whenever 
constant speed search is to be simulated, the times at sprint 
and drift speeds should be set to some large number (500 hours 
in this case) to preclude unnecessary program steps. 

The target speed is uniformly distributed between 5 
and 10 knots, with changes in course and speed taking place 
on the average of twice per hour. The target FOM is 84 db in 
the base case. 

5. Simulation Run Parameters 

Two thousand replications were run, and a maximum 
search time was set at 720 hours (30 days). ‘The random 
number seed can be any integer, and for these runs was 
selected as the time of program execution. 

6. Results and Conclusions 

Numerical results of the base case (Case A-I) and 
modified scenarios (Cases A-II, A-III, and A-IV) are shown 
in Table II. The modified scenarios reflect the following 
changes to the basic case: 

1. Case A-II: The target platform gets a 3 db gain in 
sensor performance, which results in the target and 
searcher having equal FOM at 87 db. 

2. Case A-III: The target platform reduces its radiated 
noise by 3 db which results in the target and searcher 
having equal FOM at 84 db. 

3. Case A-IV: For comparison purposes, the target 1s 
assumed to achieve both the improvements in the above 


cases, resulting in a target FOM of 87 db, anda 
Searcher FOM of 84 db. 
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The simulation results indicate that Case A-III is 
Slightly more favorable to the target than Case A-II in the 
following respects: 


1. The cumulative probability of (secure) detection by 
the searcher is smaller. 


2. A larger proportion of the searcher detections are 
direct path, resulting in decreased mean detection 
range, and increased probability of counter-detection. 


3. Less detections take place in the most distant CZ. 


4. The change in the exchange ratio in A-III (92.2) is 
greater than that for A-II (83.2). 


5. The mean time to detection is longer. 
6. The mean and median detection ranges are shorter. 

The conclusion is that a greater expected increase 
in performance will be realized if the 3 db gain is obtained 
by reducing radiated noise. Because of the better perfor- 
mance of the searcher sonar, the average slope of the 
propagation loss curve (change in prop. loss/change in 
detection range) will be less than that for the target (l.e. 
a "flatter" curve). Thus, the change in detection range 
per 3 db change would be greater for the searcher than the 
target. Therefore, reducing the searcher effectiveness (by 
quieting the target) has a larger effect on detection 
performance than increasing the target effectiveness (by 
improving the sensor) has on counter-detection performance. 

For comparison purposes, a complete set of graphical 
output from PASS for Case A-I is shown in Figure 3.4 through 


wee ane tom case A-ITIT in Figures 3.9 through 3.1/7. 
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TABLE II 


NUMERICAL RESULTS FOR CASE A EXAMPLES 


Case A-I Case A-II Case A-III Case A-IV 
PD 0.9725 0.9375 0.9210 0.8490 
PDDP 0.1496 O25 0.1900 0.1914 
PDCZ 0.8504 0.8485 0.8100 0.8086 
PDCZ1 Oa 72 0.0960 0.1482 0.0960 
PDCZ2 O. 7a22 0.7525 0.6618 0.7126 
MOES 0.9730 0.9502 0.9370 0.8732 
PCD 920015 OFOuFES O20TS0 0.0540 
PCDDP ©.3333 0.0000 0.0000 0.0093 
PCDCZ 0.6667 1.0000 1.0000 0.9907 
PEREZ 1 0.6667 1.0000 1.0000 0.9907 
PCDECZ2 0.0000 0.0000 0.0000 0.0000 
MOET 0.0080 0.0243 020310 0.0783 
ER 122). 4 39.19 30.20 aes 
(T)s Boe.) 35.2 41.8 42.9 
fr) t 26.4 17.8 32.3 33.9 
(R)s 68.9 69.7 64.6 Bue 
(R-50)s 75.6 75.6 7 oe cael 
cet 33.0 38.0 37.9 38.3 
(R=50)t a. 75 39.1 3929 8922 
Notes: 1. Due to the low number of counter-detections in Cases A-I 


and A-II (3 and 23 respectively) the variance of the target 
parameter estimates will be large. Therefore, comparison 
should be based on searcher parameter estimates. 


2. See Appendix C for explanation of abbreviations. 


56 








73.0 70.0 


PROPAGATION LOSS (0B/UPA/YD) 


120.0118.0 110.0 108.0 100.0 $8.0 0.0 6.0 


Figure 3.1. 


0.0 


SEARCHER PROPAGATION LOSS (MODEL) 


be 


10.0 2.9 W.0 @.0 SW.0 O.9 W.0 8.0 9.5 100.0 
RANGE (NAM) 


Searcher Propagation Loss Curve for Case A 
Examples. 
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Figure 3.2. Target Propagation Loss Curve for Case A 
Examples. 
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Figure 3.3. Searcher Track for Case A Examples. 
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Figure 3.4. Typical Figures-of-Merit for A-I. 
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eeoume 3.5. Cumulative Probability of Detection for A-I. 
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Figure 3.6. Progability of Detection VS. Range for A-lI. 
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Figure 3.9. Typical Figure-of-Merit for A-III. 
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Cumulative Probability of Detection for A-III. 
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Perce S212. Probability of detection Vs. Range for A-III. 
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Figure 3.17. Searcher Position at Time of Counter-Detection 
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GC. BARBER SEARCH: SOME FACTORS AFFECTING PENETRATION 
SUCCESS 


Hig Scenario Description 
A patrolling submarine (searcher) has set up a 
stationary barrier to intercept transiting submarines 
(target). Frequently submarines exhibit aspect dependent 
radiated noise patterns which may become a significant 
factor in a barrier scenario where the searcher is most 
frequently presented a bow-on aspect. It is desired to gain 
an understanding of the magnitude of this problem, and to 
investigate to what extent the distribution of targets 
across the barrier will affect intercept performance. 
#. PASS Input Data 
Table III shows the input data for the base case 
et) . 
3. Environment and Search Plan 
Figures 3.18 and 3.19 show the propagation loss 
curves for the searcher and target. Since there are no 
convergence zones, the MSEL model is not used. 
Figure 3.20 shows the search area to be a 90 nm by 
200 nm rectangle with the searcher conducting a stationary 
(back-and-forth) barrier perpendicular to the expected 
direction of target motion. The placement of the barrier 


was based on an expected median detection range of 25 nm. 
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4. Platform Parameters 
For the base case (B-I): 
1. The searcher sveed 1S constant at 10 knots. 


ey "Mme searcher FOR is 37 db and mw not target aspect 
dependent. 


3. We ainwrial target distribution across the barrier is 
wha: fomm. 


ee Sune target starts at the top of the barrier (y=200) 
(1.e. the searcher is not time late on the barrier). 


5. The maximum target course variation is plus-or-minus 
30 degrees. 


6. The target speed 1S uniform between 8 and 20 knots. 


7. The target changes course and speed (independently) 
on the average of twice per hour. 


G. Whe target FOM is 84 db. 
5. Results and Conclusions 
Numerical results for the base case (B-I) and varia- 
tions of the base case (Case B-II and B-III) are shown in 
Table IV. The variations in the base case were: 


1. Case B-II: The distribution of the initial target 
position across the ton of the barrier was altered so 
that 50% of the targets were uniformly distributed 
between 0 and 30 miles from the left boundary, and 
50% of the targets were uniformly distributed between 
09 and 30 miles from the right boundary. 


2. Case B-III: The target distribution is as in B-IlI. 
Additionally, the target was assumed to have an 
aspect dependent source level which induced the 
target aspect dependent searcher FOM as shown in 
Table V. The significance of this radiated noise 
pattern is the bow-null presented the searcher during 
barrier approach and penetration. 
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The simulation results indicate that: 
1. If the searcher is conducting a constant speed back- 
and-forth barrier, an increase in the probability of 


penetration can be realized by transiting the edge 
of the barrier. 


2. When calculating search effectiveness in a barrier 
Scenario, the aspect dependency of target source 
levels may be a significant factor. 

A complete set of graphical output from PASS for 
@eee B=iI is shown in Figures 3.21 through 3.25; for Case 


ei) im Figures 3.26 through 3.30; for Case B=III in 


Peoures 3.31 through 3.39. 
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TABLE IV 


NUMERICAL RESULTS FOR CASE B EXAMPLES 


Case B-I Case B-II Case B-III 

PD OA oss 0.6040 0.4410 

PDDP 1.0000 1.0000 lweO OOO 

MOES OreS21 0.8014 0.7152 

PCD 0.0000 0.0005 0200975 

PCDDP 0.0000 0.0000 1.0000 

MOET 0.2926 0.3949 0.5545 
ER 2-91. 2.03 1.29 

(T)s i223 12.39 12.80 
fry te none none lida 73 
(R)s 23.46 24225 PR Ae 
(R-50)s__ 22.99 24,04 Zio > 
(R)t none none 5283 
(R-50)t none none 4.69 

Notes: 1. See Appendix C for explanation of abbreviations. 


TABLE V 


SEARCHER FOM VS. TARGET ASPECT 


Target Aspect Searcher FOM 


O00 82 
O45 87 
O90 85 
135 87 
180 82 
225 87 
270 85 
aS oy 
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SERRCHER PROPAGATION LOSS (MODEL) 


73.0 70.0 





PROPAGATION LOSS (08/UPA/YD) 
90.0 85.0 


100.0 


0.9 $.0 10.0 15.0 2.0 3.0 30.0 %.0 
RANGE (NIT) 


Figure 3.18. Searcher Propagation Loss for Case B Examples. 
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TARGET PROPAGATION LOSS (MODEL) 


70.0 
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PROPAGATION LOSS (08/UPA/YD) 





0.0 5.0 10.0 18.0 20.0 25.0 w.0 
RANGE (NAM) 
Figure 3.19. Target Propagation Loss for Case B Examples. 
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SERRCHER TRACK 


Y-POSITJON (NM) 
3.0 S00 73.0 100.0 12.0 1590 175.0 200.0 


059 10.0 2.0 WO 9809 ©.9 O90 M0 W090 wW.9 
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Figure 3.20. Searcher Track for Case B Examples. 


8] 





TYPICAL F IGURE-OF -MERTI] 
MODEL: LAMBOA-SIGMA JUMP 
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OOTTED LINE: TARGET 
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Pegure 3421. Typical FPigure-of-Meérit for B-I. 
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CUMULATIVE PROBABILITY OF 
DETECTION VS. TIME 
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Perouse 3.22. Cumulative Probability of Detection for B-I. 
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CUMULATIVE PROBABILITY OF 
DETECTION VS. RANGE 
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Figure 3.23. Probability of Detection Vs. Range for B-I. 
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SEARCHER POSITION WHEN 
SEARCHER DETECTS TARGET 
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Meogure 3.24. Searcher Position at Time of Detection for B-I. 
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meaiee S625. Target Position at Time of Detection for B-l. 


86 














100.0 


FIGURE-OF -MERIT (06/UPA/YD) 
0.0 6.0 9.0 6.0 


73.0 


0.0 6.$ 


Figuee 3.26. 
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Typical Fegure=of-Merit fier BeII. 
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CUMULATIVE PROBABILITY OF 
DETECTION VS. TIME 
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Peeee 3.2/7. Cimiiative Prebabilmtyroft Deteetion for BeII. 
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CUMULATIVE PROBABILITY OF 
DETECTION VS. RANGE 
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Pere 3.28. Probability of Detection Vs. Range for 3B-II. 
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Searcher Position at Time of Detection for B-II. 


90 





TARGET POSITION WHEN 
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mere 3290. Target Position at Time of Detection for B-II. 
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TYPICAL F IGURE-OF -MERIT 
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Prquice 3.28a6 Typical Figure-of-Memit fer B-III. 


2 











CUMULATIVE PROBABILITY OF 
DETECTION VS. TIME 


0.6% 0,0 £«0.4 


0.18 0.20 


CUMULATIVE PD (UNCONDITIONAL) 
0.25 


© 
° 

"ea 
he 
- 


$.0 7.8 10.0 12.5 1$.0 17.5 
TIME (HR) 


Figure 3.32. Cumulative Probability of Detection for B-III. 
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CUMULATIVE PROBABILITY OF 
COUNTER-DETECTION VS. TIME 
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Figure 3.33. Cumulative Probability of Counter-Detection 
Com BE lil. 
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CUMULATIVE PROBABILITY OF 
DETECTION VS. RANGE 
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Miodmee 3.354. Probability of Detection Vs. Range for B-III. 
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Searcher Position at Time of Detection for 
B-ILIl. 
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Figure 3.36. 
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Target Position at Time of Detection for B-III. 
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CUMULATIVE PROBABILITY OF 
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Fegure 3.37. Pre@bability of Counter=Dection Vs. Range for 
Baer. 
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TARGET POSITION WHEN 
TARGET DETECTS SEARCHER 
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Figure 3.39. Searcher Position at Time of Counter-Detection 
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D. DEFINITE RANGE LAW APPROXIMATIONS 


l. Area Search: Best Patrol speed to Evade a SY Sitematic 


Searcher 
a. Scenario Description 

For Case C-I the searcher is given a definite 
range law sensor (detection range of 20 nm), and the 
stationary target (uniformly distributed over the search 
area, A) has no capability to counter-detect. The searcher 
SyStematically sweeps the search area at a constant speed 
(v) with track spacing equal to the detection range. See 
Figure 3.24. After 88 hours of search (at 5 knots) the 
searcher has come within one detection range of every point 
in the search area, or has completely swept the area. 

In Case-II target motion iS introduced. The 
target travels at a constant speed of 5 knots, and changes 
course randomly on the average of twice an hour. 

me Paes Imput Data 

The input data for Case C-I is shown in Table VI. 
Note that in these examples figure-of-merit is chosen to 
give detection ranges based on propagation loss curves shown 
in Figures 3.40 and 3.41. 

c. Results and Conclusions 

Numerical results for Case C-I and Case C-II are 
Mom in Table VII. Graphical results for Case C-I are 
Siem in Figures 3.43 through 3.46, and for Case C-II in 


Figures 3.47 through 3.49. 
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A comparison of the results shows: 


1. The times to detection are roughly uniform for Case 
C-I as can be seen in the very nearly constant slope 
of Figure 3.43. Note the departure from exhaustive 
search results (manifested by the jump at time = 0 in 
Figure 3.43), which is due to the significant prob- 
ability of detection at the instant the search 
Starts. Also, "bumps" in the cumulative probability 
of detection curve occur at turn points, and repre- 
sent the combined effects of the searcher having a 
Slgnificant portion of his detection circle outside 
the search area, and sweeping area previously swept, 
thus resulting in a low detection rate. 


2. Introducing target motion in Case C-II skews the 
distribution of time to detection to the right (see 
Figure 3.47), resulting in an increased mean time 
to detection. 





Conclusion 2 above raises the question: What is 
the "best" speed to evade? The results of Case C-II indicate 
that the best speed 1S not zero, if the criteria 1S mean 
time to detection. Before proceeding, the following terms 
are defined: 

1. Speed Ratio: The ratio of target speed (u) to 
searcher speed (v). Note that if the speed ratio is 
zero, then the target 1s stationary. 

2. Mean Time Ratio: The ratio of the mean time to detec- 
tion (for a given speed ratio) to the mean time of 
detection when the speed ratio 1S zero. 

3. Median Time Ratio: The ratio of the median time to 
detection (for a given speed ratio) to the median 
time to detection when the speed ratio 1s zero. 

To investigate the problem of best evasion 
speed, PASS was run for various search speeds ( v = 5, 10, 


15, 20 knots) and various speed ratios at each search speed, 


under case C conditions (definite range law sensor). 


Moz 














Figures 3.50 and 3.51 show the results of these 
Simulations. Note that whenever the time ratio is greater 
than 1.0, the time to detection 1S greater than the time to 
detection for a stationary target. Figure 3.50 shows the 
mean time to detection peaks at speed ratios between 1.0 and 
1.4, depending on the searcher speed. Figure 3.51 shows the 
median time to detection is relatively insensitive to speed 
ratios below 1.0, and drops off rather sharply for speed 
ratios above 1.0. The conclusion drawn is that the best 
Speed to evade the searcher is about equal to the searcher 
Speed. If the searcher speed is not known, the best tactic 
would be to limit maximum patrol speed to the most likely 
minimum searcher speed. It is important to note that the 
best patrol speed is not necessarily the minimum achievable 


speed. 


103 











0°O = (€)WWSIS 0-0 = (c)moIs 0°O = (T)WWOIS TOO°O = (€)VaanvT 

TOO°O = (c)Vadnvl TOO°O = (T)VYaanvl 00°88 = XWNL O00S = dayn G€cc6 = Gags 

OO°Z = DOLY TOO°O = OSLY 0°O = XWNLS 0°O = NIWLS 00°OL SLWOA 00°OL GINO 
00°OOT = SL 00°O00T = GL 00°S = SOS 00°S = GOS 00° 6 SOWNOA 00°P6 = GOWOI 


>SToVoweAeCG uNY pue WIAOCTRAe Td HutTutTewsy 


00°OST = (OT) TL 00°S6 = (OT) LY 00°O0¢T = (OT)TO 00°S6 = (O0T)OU 
0O0°9PT = (6)TL 00°S8 = (6) Lu OO°9TT = (6)T0 00°S8 = (6)O" 
OO°9ET = (8) TL 00°09 = (8) Lu 00°s0T = (8)7T0 00°09 = (8)O¥ 
OO°OET = (L) TL 00°St = (L) Lu OO°PpOT = (ZL)TO 00°SP = (L)OU 
OO°STT = (9)TiL oo°O€ = (9) LY 00°96 = (9)T0 OO°O€ = (9)O0' 
00°60T = (S)TL 00°O? = (S) Lu 00°P6 = (S)TO 00°02 = (S)O"U 
O00°SOT = (b) IL OO°OT = (fF) LY OT°€6 = (%)TO OO°OT = (7)O' 
OO°POT = (€) TL OS°L = (€) LY 06°c6 = (€)TO OS°L = (€)0OU 
00°S6 = (C)IL 0O0O°S = (2) LY 00°S8 = (2)TO 00°S = (2)O"U 
0O0°OL = (T) TL OS°c = (1) Le 00°OL = (T)TO OS°e = (T)OU 
>SsSoTJT uoTtjzebedorg 3ehrAeL :SsoyT uoTtjzebedorg AsYyOAe|S 
OO°O¢T = (9)dAzA 00°02 = (9)dx 
0°O = (S)dA 0O0°O0c = (S)dX 
0°O = (P)dA 00°09 = (%)dxX 
OO°O7T = (€)dA 00°09 = (€)dx 
OO°O?T = (2)dA OO°OOT = (2) dX 
0°O = (T)d&A o0°OOT = (T)dxX O00°O@T = XWWA 
1 = @ 9 = dN 00°O¢@T = XWWX 
>SQuUTOg AoYyOUY YOeAZL AeyoAess >SuUOTSUSWTG ePoeAY YOAeRDS 


I-O FSWO HOH VLVd LNdNI 


TA AIL 





104 





TABLE VIT 


NU@MERLECAL RESULTS FOR CASES C-I AND C=-II 


PD 
EDDP 
MOES 


Pep 
PEDDP 
MOET 


ER 


CLs 
Vartr)s 
T (min) 


T (max) 


Notes: 


Case C-I 


1.0000 
1.0000 
1.0000 


0.0000 
0.0000 
0.0000 


infinite 


40.12 
TOG ay 
GROG 
87.24 


Case C-II 


ORg292 
1.0000 
0.9994 


0.0000 
Q0.0000 
0.0002 


4496.00 


43.84 
Jim. 38 
0.00 
302.07 


1. See Appendix C for explanation of 


abbreviations. 
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SEARCHER PROPAGATION LOSS (MODEL) 


3.0 70.0 


PROPAGATION LOSS (D8/UPA/Y0) 


120.011$.0 $10.0 108.0 100.0 95.0 9.0 05.0 0.0 


6.0 10.0 0.0 NWO 0.0 0.2 8.8 0.0 0.0 8.8 100.0 
RANGE (NEM) 


Figure 3.40. Searcher Propagation Loss for Case C Examples. 
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TARGET PROPAGATION LOSS (HODEL} 
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Figure 3.41. Target Propagation Loss for Case C Examples. 
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Figure 3.42. Searcher Track for Cases C-I and C-II. 
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CUMULATIVE PROBABILITY OF 
DETECTION VS. TIME 
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Pro@se 3443. Cumallative Probability of Detection for C-lI. 
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TARGET POSITION WHEN 
SEARCHER DETECTS TARGET 
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Figure 3.44. Target Position at Time of Detection for C-I. 
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Figure 3.45. Searcher Position at Time of Detection for C-I. 
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CUMULATIVE PROBABILITY OF 
DETECTION VS. RANGE 
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Figure 3.46. Probability of Detection Vs. Range for C-I. 
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Figure 3.47. Cumulative Probability of Detection for C-II. 
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SEARCHER POSITION WHEN 
SERRCHER DETECTS TARGET 
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Figure 3.48. Searcher Position at Time of Detection for 
Caan, 
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MEAN TIME RATIO VS. SPEED RATIO 
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MEDIAN TIME RATIO VS. SPEEO RATIO 
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E. CONCLUSIONS AND RECOMMENDATIONS FOR FURTHER WORK 

In the preceeding examples, it was intended to demon- 
strate the capabilities and flexibility of the PASS program. 
Of significance is the fact that an experienced operator, 
using the interactive data entry procedures, could run all 
of the example problems in about four hours (including time 
for hard-copy printing and graphics production). These runs 
represent 274,000 individual simulation replications. 

The utility of the graphical output should be evident. 
Frequently the pictorial displays provide insight that is 
obscured by a large quantity of numbers. Additionally, the 
very shape of certain graphical results inspire questions 
the answer to which are critical to a thorough analysis. 

Improvements or additions to PASS might include: 

1. Multiple sensor capability by both platforms. 


2. Multiple search platforms (e.g. searcher, target, 
and consort). 


3. Post detection analysis (e.g. break contact, 
approach and attack, target motion analysis). 


4. Gaming capability (real time control over platform 
Mei Lon) . 


5. Weapon employment capability. 
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APPENDIX A 


PASS USERS GUIDE 


me. INTRODUCTION 

This User's Guide is a set of instructions designed to 
allow the execution of program PASS on the NPGS computer by 
users unfamiliar with the program code. The user should be 
familiar with the concepts of the simulation model as 
presented in Chapter I, and have available for reference the 
variable list in Appendix B. 

1. Data File Management 

When PASS is run, an input data file may, or may not 
exist. If no data file exists, the easiest way to create one 
is to execute program PASS and input data during the inter- 
active session. If a data file exists it must have the 
filename PASS and filetype DATA. If a data file is inter- 
actively created, it will be in file PASS DATA at the 
completion of the run of PASS. 

At the start of each run, if PASS DATA exists, i1t 1S 
copied into file PASS HOLD. If changes are made to the input 
data, file PASS DATA will reflect these changes, and PASS 
HOLD will contain the unchanged data. If no changes were 
made, PASS HOLD and PASS DATA will be identical. To save the 
data in PASS HOLD, it must be renamed before the next run of 
PASS. PASS DATA will remain unchanged as long as no changes 


are made to the data during the interactive session. 
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Upon completion of each PASS run a file PASS OUTPUT 
is created. Depending on the option selected by the user, 
this file will contain some, or allofthe following: 


l. A complete record of the input data, or a user 
Specified run identification number. 


2. A summary of the results of the simulation. 

3. Sectioning results of times and ranges of detection 
and counter-detection (requires a minimum of 500 
data points). 

4. Histogram and statistics of times and ranges of 
detection and counter-detection (requires a minumum 
of 10 data points). | 

If this file is to be saved, it must be renamed or printed 
before the next run of PASS. 
Upon completion of each run of PASS a file PASPLT 
DATA is created. This file provides the data, in the proper 
format, to be used by the FORTRAN program PASPLT. Program 
PASPLT uses the DISSPLA graphics — resident on the NPGS 
computer to produce a graphical display of the simulation 
input parameters and results. 
a Required Input Data 
PASS requires the following data, which the user 
should have ready for interactive data input if a complete 
PASS DATA file does not exist: 
im Saze of the search area rectangle. 
2. Searcher track anchor points and return point. 


3. Searcher propagation loss data (direct path and CZ'‘s). 


4. Target propagation loss data (direct path and CZ's). 
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7. Seeeet Gaatial X—dilsitributzen if met umiform (barrier 
scenario only). 


G. Searcher figure-of-merit at sprint and drift speeds 
tera “function of target aspect,.if target source 
level is aspect dependent). 


7. Searcher sprint and drift speeds. 


8. Searcher sprint and drift times. 
aelenget fvgure-of=merit against sprinting and drifting 
Searcher. 


10. Target speed range. 

ll. Rate of target zigs (course and speed changes). 
12. Acoustic fluctuation rate parameters. 

13. Acoustic fluctuation scale parameters. 


14. Run parameters (random number seed, number of repli- 
cations, maximum search time for each replication). 


During the interactive input of data, the program 
does some error checking. However, it is far from impossible 
to enter bad data, or incorrectly formatted data. If the 
program fails to run as expected, the most probable cause is 
bad data. 

cewenie Execubive Program 

The executive program PASS EXEC is necessary to run 
the FORTRAN program PASS. The executive program conducts 
data file management and starts PASS. Prior to invoking the 
executive program, PASS must be compiled in either FORTRAN-G 
Or FORTRAN-H with the text file on the A-DISK. The 
FORTRAN-H compiler is recommended as it results in a faster 


program run time. 
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B. PROGRAM PASS EXECUTION 

Each of the following sections describes a terminal 
option. Each option presents a new menu on the terminal 
screen. Depending on your response to each option menu, you 
will proceed to the next option, or you will be prompted to 
enter data or make some decision regarding the nature of the 


Simulation. 


im Starting the Program: Invoking the Executive 
Program 


MS Stare program PASS, ENTER: PASS 
The screen will contain a quick review of data-file manage- 
ment and the first OPTION. 
2. OPTION NO. 1: Terminate or Proceed 
This option allows you to "gracefully" terminate 
the program. 
To terminate the program, ENTER: l 
To proceed with execution, ENTER: 2 
If you terminate now, and you had an input data file PASS 
DATA, it will have been copied into file PASS HOLD, and PASS 
DATA will be empty. 
foe fect ER OPTION: Accept Program Defaults 
This option allows you to bypass the interactive 
sections and proceed directly to the simulation. Bypassing 
requires a complete and properly formatted input data file, 


and the acceptance of the program defaults. 
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To list the program defaults, ENTER: 1 
To accept the program defaults, ENTER: 2 
To commence interactive options, ENTER: 3 


4 OPTION NO. 2A: Signal Integration Model 


This option allows you to select a threshold 


crossing model without signal integration, or the MSEL 


integration model. It 1s recommended that the MSEL model 


be used if strong convergence zones are present. 
To GikeGtE no signal immtegtation, ENTER: 1 
To select the MSEL integration model, ENTER: 
5. OPTION NO. 2: Acoustic Fluctuation Model 
This option allows you to select either a 
Lambda-Sigma Jump or a Gauss-Markov error process to 


acoustic fluctuations. 


To select the Lambda-Sigma Jump model, ENTER: 


To select the Gauss-Markov model, ENTER: 2 
6. OPTION NO. 3A: Existence of File PASS DATA 
The program must be "told" whether or not to 
data from an existing data file. 
If PASS DATA is not on your A-DISK, ENTER: l 
tt PASS DATA ws on your A=DISK, ENTER: 2 
7. OPTION NO. 3: Search Area Size 


Variables: XMAX, YMAX 


2 


mode 1 


read 


This option allows you to specify the size of the 


search area in which the target is confined. The area dimen- 


sions are in nautical miles. 
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To review existing data, ENTER: 1 

To accept existing data, ENTER: 2 

To change existing data, ENTER: 3 

If (1) 1s entered, the existing data (in file PASS 
DATA on your A-DISK) will be presented on the screen for 
BeView. To return to the option menu to accept or modify 
this data, ENTER: 0. 

If (2) 1s entered, the existing data will be read 
into PASS, and the next option menu will be presented on 
the screen. 

If (3) 1s entered, prompting messages will be on the 
terminal screen to allow you to enter new data. When data 


entry 1s completed, you must return to the option menu (to 


meview, accept, or further modify data). To return, 
BMPR: O. 
8. OPTION NO. 4: Searcher Track 


Variables: NP, KP, XP(I), YP(I) 

This option allows you to enter the searcher track 
anchor points (2 to 50 points), and the return point. The 
dimensions of the anchor points are nautical miles from 
or 0) . 

To review existing data, ENTER: l 

To accept existing data, ENTER: 2 

To change existing data, ENTER: 3 

If (1) is entered, the existing data (in file PASS 


Pere On your A-DISK) will be presented on the screen 
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for review. To return to the option menu to accept or 
modify this data, ENTER: 0. 

If (2) 1s entered, the existing data will be read 
into PASS, and the next option menu will be presented on 
the screen. 

If (3) is entered, prompting messages will be on the 
terminal screen to allow you to enter new data. When data 
entry is completed, you must return to the option menu (to 
beview, accept, or futher modify data). To return, 

ENTER: O. 
Se (CeTMIONTNO...5:  Séeaecher Direct Path Propagation Loss 

Variables: MO, RO(I), OL(I) 

This option allows you to input the propagation loss 
of the searcher as a function of range. The number-of points 
should be between 2 and 20, enSuring that there are suffi- 
cient points to adequately describe the propagation loss 
Curve. Range is in nautical miles, and propagation loss in 
decibels (db). 

To review existing data, ENTER: l 

To accept existing data, ENTER: 2 

Ho change exesting datayeENTER: 3 

If (1) is entered, the existing data (in file PASS 
DATA on your A-DISK) will be presented on the screen for 
review. To return to the option menu to accept or modify 


Eames Gata, ENTER: O. 
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If (2) 1s entered, the existing data will be read 
into PASS, and the next option menu will be presented on 
the screen. 

If (3) 1s entered, prompting messages will be on the 
‘terminal screen to allow you to enter new data. When data 
entry 1s completed, you must return to the option menu (to 
review, accept, or further modify data). To return, 

GWRER: 0. 
HO> OPTION NO. 6: Target Direct Path Propagation Logs 

Variables: MT, RT(I), OL(I) 

Theg Option 1s samalar to OPTION 5, but here you 
enter the propagation loss data for the target. 

To review existing data, ENTER: l 

To aceept exasting data, BNPER: 2 

To change existing data, ENTER: 3 

If (1) is entered, the existing data (in file PASS 
DATA on your A-DISK) will be presented on the screen for 
review. To return to the option menu to accept or modify 
fears Gata, ENTER: OQ. 

If (2) is entered, the existing data will be read 
into PASS, and the next option menu will be presented on 
the screen. 

If (3) is entered, prompting messages will be on the 
terminal screen to allow you to enter new data. When data 
entry is completed, you must return to the option menu (to 


mevnew, accept, or further modify data). To return, ENTER: QO. 
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[3. OPPICGH NO. 7: Searcher Convergence Zones 
Vana@bles: MGZO, ROZO(I,J), GZLO(L) 


This option allows you to input the number of 
Searcher convergence zones (0 to 5), the width of each CZ, 
and the propagation loss in each CZ. The CZs are modeled as 
inverted square-wells superimposed on the DP propagation 
loss curve. You will have to make a subjective interpreta- 
tion of the actual CZ data to determine the range to the 
inner and outer rings of the square-well, and the effective 
propagation loss therein. Ranges are in nautical miles, 
amd propagation lioss in db. 

To review existing data, ENTER: l 

To accept existing data, ENTER: 2 

Techamge existing data, ENTER: 3 

If (1) is entered, the existing data (in file PASS 
DATA on your A-DISK) will be presented on the screen for 
review. To return to the option menu to accept or modify 
Paes Gata, ENTER: QO. 

If (2) is entered, the existing data will be read 
into PASS, and the next option menu will be presented on 
the screen. 

If (3) is entered, prompting messages will be on 
the terminal screen to allow you to enter new data. When 
data entry is completed, you must return to the option menu 
[eo review, accept, or further modify data). To return, 


ENTER: QO. 
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ha: OPTION NO. .8: Target Convergence Zones 
Variables: NG2T, REZT(1,73);, CGZBT (1) 


Thag Option bs similar to OPTION 7, but hexe you 
enter the data for target convergence zones. 

To review existing data, ENTER: l 

To accept existing data, ENTER: 2 

To change existing data, ENTER: 3 

If (1) is entered, the existing data (in file PASS 
DATA on your A-DISK) will be presented on the screen for 
review. To return to the option menu to accept or modify 
Cmes Gata, ENTER: 0. 

If (2) is entered, the existing data will be read 
into PASS, and the next option menu will be presented on 
the screen. 

If (3) is entered, prompting messages will be on 
the terminal screen to allow you to enter new data. When 
data entry is completed, you must return to the option 
menu (to review, accept, or further modify data). To 
return, ENTER: O. 

13. OPTION NO. 9: Area or Barrier Search 

This option allows you to select either an area or 
barrier search scenario. If area search is selected, the 
next option will be 1LOA. If barrier search is selected, the 
next option will be 9A. 

To select area search scenario, ENTER: l 


To select barrier search scenario, ENTER: 2 
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4. OPTION NO. 9A: Initial Terget Lateral Distribution 
Variables: NBINS, XBIN(I), PBIN(TI) 


This option allows you to specify the lateral 
@@etrtbution of target starting position (barrier scenario 
only). Distances are in neil anit miles from X = Q. 

To distribute the initial positions uniformly on 
(O7oaAxX), ENTER: 1. 


To specity the lateral distribution, EBNIER: 2 


mo CrmION NO. 1O0A: Searcher FOM Dependent on, Target 
Aspect 


Variables: NLS, BRG(I), FOMBD(I), FOMBS(TI) 


This option allows you to specify the searcher 





| figure-of-merit (at drift and sprint speeds) as a function 

| of the relative bearing from target to searcher. This 1s 

| equivalent to a target source level dependent upon aspect 
angle. Angles are in degrees, and FOM in db. 

To specify searcher FOM not target aspect depen- 
Gent, ENTER: 1 

To specify searcher FOM target aspect dependent, 
Pee R: 2 

If (1) is selected, OPTION 10 will be presented next. 

Tf (2) is selected, after data entry, OPTION 11 will be 


presented next. 
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mS. OPTION NO. 10: Searcher Figure-of-Merit 
Variables: FOMOS, FOMOD 


This option allows you to specify the searcher 
fFigure-of-merit against the target when the searcher is at 
Sprint or drift speeds. FOM IN db. 

TO review existing data, ENTER: 1 

To accept existing data, ENTER: 2 

To change existing data, ENTER: 3 

If (1) is entered, the existing data (in file PASS 
DATA on your A-DISK) will be presented on the screen for 
review. To return to the Option menu to accept or modify 
this data, ENTER: 0. 

If (2) is entered, the existing data will be read 
into PASS, and the next option menu will be presented on 
the screen. 

If (3) is entered, prompting messages will be on 
the terminal screen to allow you to enter new data. When 
data entry is completed, you must return to the option 
menu (to review, accept, or further modify data). To 
Boeemn, ENTER: Q. 

OL eee NO. llcenSeancher  Gpuunt. and. Drift Speeds 

Variables: SOS, SOD 

This option allows you to select the search speed 
(Sprint and drift). For constant speed search, set SOS = 


SOD. Speeds in knots. 
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To review existing data, ENTER: 1l 

To accept existing data, ENTER: 2 

To change existing data, ENTER: 3 

Tf (1) is entered, the existing data (in file PASS 
DATA on your A-DISK) will be presented on the screen for 
review. To return to the option menu to accept or modify 
thas data, ENTER: O. 

If (2) is’entered, the existing menu will be read 
into PASS, and the next option menu will be presented on 
the screen. 

If (3) is entered, prompting messages will be on 
the terminal screen to allow you to enter new data. When 
data entry is completed, you must return to the option 
menu (to review, accept, or further modify data). To 
return, ENTER: 0O. 

mee |6OP IEC? MO. 12: Gearcher Sprint and Drift Times 

Variables: TS, TD 

This option allows you to specify the time the 
Searcher spends at sprint and drift speeds. Time in hours. 

To review existing data, ENTER: l 

To accept existing data, ENTER: 2 

To change existing data, ENTER: 3 

If (1) is entered, the existing data (in file PASS 
Mets On your A=DISK) will be presented on the screen for 
review. To return to the option menu to accept or modify 


this data, ENTER: 0. 
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If (2) is entered, the existing data will be read 
into PASS, and the next option menu will be presented on the 
screen. 

If (3) is entered, prompting messages will be on the 
terminal screen to allow you to enter new data. ten Jews 
entry is completed, you must return to the option menu (to 
review, accept, or further modify data). To return, 

ENTBR: O. 
Jae OPREONTNO. 13: Target Figure-of-Merit 

Variables: FOMTS, FOMTD 

This option allows you to specify the target 
Sigure-of-merit against the searcher when the searcher is at 
Sprint or drift speed. FOM in db. 

To review existing data, ENTER: 1 

To accept existing data, ENTER: 2 

To change existing data, ENTER: 3 

If (1) is entered, the existing data (in file PASS 
DATA on your A-DISK) will be presented on the screen for 
review. To return to the option menu to accept or modify 
tas data, ENTER: QO. 

If (2) is entered, the existing data will be read 
into PASS, and the next option menu will be presented on the 
screen. 

If (3) is entered, prompting messages will be on 
the terminal screen to allow you to enter new data. When 


data entry is completed, you must return to the option 
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imam (to review, accept, or further modify data). To 
Peatuaem, ENTER; 0. 
20. OPTION NO._14: Target Speed Range 

Variables: STMIN, STMAX 

This Option allows you to specify the range of 
target spped. Target speed is distributed uniform on 
[uy oTMAX). Speed in knots. 

To review existing data, ENTER: l 

To accept existing data, ENTER: 2 

To change existing data, ENTER: 3 

If (1) is entered, the existing data (in file PASS 
DATA on your A-DISK) will be presented on the screen for 
review. To return to the option menu to accept or modify 
eaas$ data, ENTER: QO. 

If (2) is entered, the existing data will be read 
into PASS, and the next option menu will be presented on 
the screen. 

If (3) is entered, prompting messages will be on 
the terminal screen to allow you to enter new data. When 
data entry is completed, you must return to the option menu 
(to review, accept, or further modify data). To return, 
Bieter: 0. 

21. OPTION NO. 15: Rate of Target Zigs 
Variables: _ RTCC, RTSC 
TMrusmopelon allews you to specify the rate of target 


course and speed changes. Rates are "per-hour”". 
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To review existing data, ENTER: l 

To accept existing data, ENTER: 2 

To change existing data, ENTER: 3 

If (1) 1s entered, the existing data (in file PASS 
DATA on your A-DISK) will be presented on the screen for 
review. To return to the option menu to accept or modify 
weaSea@ata, ENTER: 0. 

If (2) 1s entered, the existing data will be read 
into PASS, and the next option menu will be presented on the 
screen. 

If (3) 1s entered, prompting messages will be on the 
terminal screen to allow you to enter new data. When data 


entry 1S completed, you must return to the option menu (to 


review, accept, or further modify data). To return, 
ERTER: O. 
22. OPTION NO. 16: Acoustic Fluctuation Rate Parameters 


Variables: ALAM(T) 

This option allows you to specify the rate param- 
eters for the acoustic fluctuation process. Rates are 
"per-hour". 

To review existing data, ENTER: 1 

To accept existing data, ENTER: 2 

To change existing data, ENTER: 3 

If (1) is entered, the existing data (in file PASS 


DATA on your A-DISK) will be presented on the screen 
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for weveew. TO retumn to the option menu to accept or modify 
mas Gata, ENTER: 0. 

If (2) is entered, the existing data will be read 
into PASS, and the next option menu will be presented on 
the screen. 

If (3) is entered, prompting messages will be on 
the terminal screen to allow you to enter new data. When 
data entry 1s completed, you must return to the option menu 
(EO neview, accept, or further modify data). To return, 
EPeeeR: 0. 


23. OPTION NO. 17: Acoustic Fluctuation Scale 
Parameters 


Variables: SIGMA(T) 

This option allows you to specify the scale param- 
eters for the acoustic fluctuation process. Scale in db. 

To review existing data, ENTER: l 

Fo aceept existing data, EN@ER: 2 

To change existing data, ENTER: 3 

Tf (1) is entered, the existing data (in file PASS 
DATA on your A-DISK) will be presented on the screen for 
review. To return to the option menu to accept or modify 
this data, ENTER: QO. 

Tf (2) is entered, the existing data will be read 
into PASS, and the next option menu will be presented on 


the screen. 
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If (3) is entered, prompting messages will be on 
the terminal screen to allow you to enter new data. When 
data entry is completed, you must return to the option menu 
wee review, accept, or futher modify data). To return, 
BN@eER: 0. 

2@5 OPTION NO. 18: Run parameters 

Variables: SEED, NREP, TMAX 

This option allows you to select parameters that 
govern the simulation execution and random number generation. 
Time in hours. 

TO review exlsting data, ENTER: 1 

ie aeceept existing data, ENTER: 2 

To change existing data, ENTER: 3 

If (1) is entered, the existing data (in file PASS 
DATA on your A-DISK) will be presented on the screen for 
review. To return to the option menu to accept or modify 
this data, ENTER: 0. 

If (2) is entered, the existing data will be read 
into PASS, and the next option menu will be presented on 
the screen. 

If (3) is entered, prompting messages will be on 
the terminal screen to allow you to enter new data. When 
data entry is completed, you must return to the option 
menu (to review, accept, or further modify data). To 


return, ENTER: 0O. 
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a=» OPTION NO. 19: ‘Input Deiea to Output File 


This option allows you to either send all of the 
maput data to the output file, or to identify the output 
file with a run identification number. 

To send all input data to the output file, 
BwIpR: 1 

TO Supress input data echo and select a run ID, 
BNTER: 2 

Zoe OPTION NO. 205% Statistical Amalysis 

This option allows you to select the amount and 
type of statistical analysis to be done on the simulation 
results. 

For complete statistical analysis (sectioning, 
histograms, and parameter point estimates, ENTER: l 

Tomeomutesect OonanggveENTER: 2 


For no statistical analysis, ENTER: 3 


C. NORMAL RUNNING INDICATION 

Immediately following OPTION 20, the screen will 
display a running total of the number of replications 
completed. This display is updated every 200 replication, 
which allows the user to follow the progress of program 
@ugeltdionm. Upon completion of the last replication, the 
Screen will be cleared, and the simulation quick-summary 


will be presented. 
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D. GRAPHICAL PRESENTATION OF RESULTS: PROGRAM PASPLT 

The FORTRAN program PASPLT is intended to be used in 
conjunction with the DISSPLA graphics system to produce 
graphical representations of the results of a run of program 
Peoo. To run program PASPLT you will need: 


? Data file PASPLT DATA on your A-DISK. PASPLT DATA is 
automatically created by running program PASS. 


2. Dual terminal graphics capability (TEKTRONIX 618 and 
thermal printer). 


3. Program PASPLT FORTRAN-G or FORTRAN-H compiled on 
youre A-DISK. 


To produce graphical output, invoke the graphics package 
by entering: DISSPLA. Two responses to the DISSPLA execu- 
tive program prompts are required: 


1. You must identify the program PASPLT as the FORTRAN 
program to be used with DISSPLA. 


2. You must define the data file as follows: FILEDEF 10 
DISK PASPLT DATA. 


The above is done in response to prompting messages produced 
on the screen after DISSPLA is entered. 

Each plot is presented on the TEKTRONIX 618 screen, and 
hard copies may be made by pressing the HARD COPY button at 
the bottom of the screen. To produce the next plot, hit 
ENTER on the keyboard. If any messages (other than error 
messages) appear on the IBM screen, clear the screen until 


the next plot appears on the TEKTRONIX screen. 
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APPENDIX B 


DESCRIPTION OF VARIABLE NAMES USED IN PASS 


The following list of major variable names used in PASS 
is arranged alphabetically. Those variable names with an 
asterik (*) indicate that an understanding of their function 
1S required to run PASS in conjunction with Appendix.A (PASS 
User's Guide). Following each variable name is the data type 
and physical dimension, if any, and then a verbal description 
of the variable use. For input data, restrictions (e.g. 
variable size, entry procedures) are listed under "RESTR:". 
Some variables used for temporary storage, loop counters, 

Or pointers are not listed as their function should be 


apparent to the reader of the code. 


PASS VARIABLE NAMES: 


1. AFL(I)/real/decibels/..... Value of the Ith acoustic 
fluctuation process. I=1,2,3. 

Zee ee eM /logical/..... If the Gauss-Markov error process 
is to be used, this variable is TRUE. Otherwise, it 
is false. 

oe Mees /logical/..... If the Lambda-Sigma Jump error 
process is used, this variable is TRUE. Otherwise, 


ite 1S FALSE. 

4. ALAM(I)/*/real/inverse hours/..... Rate parameter for 
@fe. Wea fluctuation process, I1=1,2,3, (input). 
RESTR: positive. 


wm, ABBA/real/..... Multiplier for the GMA process. 


LB) 





nO. 


al. 


2 


PS. 


14. 


IPD . 


ro. 


ANG1/real/degrees (external), radians(internal)/..... A 
limit of target courses allowed in the barrier scenario. 
In the simulation, courses are measured CCW from the 
X-axis. In the barrier scenario, the target base 
course iS in the negative Y-direction, or psuedo- 
course of 270-degrees. ANGI] = 270 - ANG2. 


ANG2/*/real/degrees (external), radians(internal)/..... 
Angular variation from base course (ANG1l) allowed for 
the target in the barrier scenario. The target courses 
are distributed uniformly over psuedo-course 270+ANG2. 
See ANG] (input). RESTR: non-negative, not greater 
than 90-degrees. 


PSPBGT/logical/..... If the target source level is 
aspect dependent, this variable is TRUE. Otherwise, 
1t is FALSE. 


BETA/real/..... Multiplier for the GMA process. 


BREL/real/radians/..... Relative bearing from the 
target to the searcher. 


BRG(1) /*/real/degrees (external), radians (internal) / 
Ms ees Relative bearing associated with target aspect 
related searcher figure-of-merit, I = l,...,NLS, 
Gumput). RESTR: BRG(1l) must be 000.0 (directly on 
the bow), and all subsequent BRG(I) must be in 
ascending order, such that BRG(I) < BRG(I+1), and 
BRG(I) max = 360.0. 


fie le? we eI aNnteger/..... Counters for STACK1,2,3,4. 


G@enm/meal/..... X-component of a unit vector in the 
direction of target motion. 


GMO (1) /*/real/decibels/..... Propagation loss in the 
Ith. searcher convergence zone (input). If the number 
of searcher convergence zones (NCZO) is zero, then 
tes data us not input. Otherwase, I=l,.«...,NCZO. 
PESPR: positive. 


CZLT(I)/*/real/decibels/..... Propagation loss in the 
Ith target convergence zone (input). If the number 
of target convergence zones (NCZT) is zero, then this 
data in not input. Otherwise, I=l1,...NCZT. 

RESTR: positive. 


Dl/ceoal/nawtrcal miles/..... Distance from current 


Searcher position (X0,YO) to the projected searcher 
pesmeton (XO, YOT) . 
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D2/real/nautical miles/..... Distance from current 
searcher position (X0,YO) to the end of the current 
search leg. 


Bier (i) /cteal/nawtacal milesi/..... The distance from 
Searcher track anchor point I to anchor point (I+1). 
The length of the Ith leg of the search pattern. 


Bids) /cead/as sss X-component of the unit vector in the 
direction of searcher motion from anchor point I to 
anenor pommt (I*+1). 


Dial) /rea@l/... Y-component of the unit vector 
pointing in the direction of searcher motion from 
g@mehOr.point I to amchor point (I+1). 


P@n@S/ logical / wo... If the user chooses to have the 
input data sent to the output file, this variable is 
Pawn. Orhemwrse, 1t 1S PALSBs 


BNOT /logu@al/s..s. If 5 hours of figure-of-merit data 
for the target has been recorded, this variable is 
TRUE. Otherwise, it 1s FALSE. 


BNUPF/logical/..... If 5 hours of figure-of-merit data 
for the searcher has been recorded, this variable is 
TRUE. Otherwise, it 1s FALSE. 


MO@(1) /real/fdecibels/..... A sample of the searcher 
figure-of-merit for the first five hours of a repli- 
cation used to provide a representative plot of FOM 
as a reference. This variable is not used computa- 
tionally in the program. 


FOMBD(I)/*/real/decibels/..... Searcher figure-of- 
merit when at drift speed and with a relative bearing 
from the target to the searcher of BRG(I), where 


io lee. ,teuo, (input). RESTR: positive. 

FOMBS (I) /*/real/decibels/..... Searcher figure-of- 
merit at sprint speed at the corresponding relative 
bearing from target to searcher of BRG(I), I =1l,..., 
NLS, (input). RESTR: positive. 
FOMO/real/decibels/..... Current searcher figure-of- 
merit. 

FOMOD/*/real/decibels/..... The searcher figure-of- 
merit at drift speed (input). RESTR: positive. 
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FOMOS /*/real/decibels/..... The searcher figure-of- 


merit at sprint speed (input). RESTR: positive. 
FOMT/real/decibels/..... Current target figure-of- 
merit. 

FOMTD/*/real/decibels/..... The target figure-of-merit 


when the searcher is at drift speed (input). 
RESTR: positive. 


FOMTS/*/real/decibels/..... The target figure-of-merit 
when searcher is at Sprint speed (input). 
RESTR: positive. 


Fae /logical/..... This variable is TRUE if the user 
Specifies a data file, PASS DATA, presently exists on 
his A-disk. Otherwise, it is FALSE. 


FT (I) /real/decibels/..... A sample of the target 
figure-of-merit for the first five hours of a repli- 
cation used to provide a plot of representative FOM 
as a reference. This variable is not used computa- 
tionally elsewhere in the program. 


GO4IT/logical/..... This? vaweilable is TRUE if tie user 
Specifies that the interactive data input is to be 
bypassed and input data is to be read directly from 
an existing file on A-disk, PASS DATA. Otherwise, 
wt is FALSE. 


B@ZO7integer/..... Identifies the number of the 
convergence zone in which the searcher CZ detection 
took place. 


nemm/integer/..... Tdentifies the number of the 
convergence zone in which the target CZ counter- 
detection took place. 


mimrec/logical/..... This varlable is TRUE if the MSEL 
model is to be used. Otherwise, it 1s FALSE. 


KBAR/integer/..... Indicator variable. KBAR = 1 if 
the area search@scenario is used. KBAR = 2 if the 
barrier search scenario is used. 


KCZO/integer/..... The number of the most distant CZ 
in which the searcher can currently make a detection. 
If KCZO = 0, no CZ detections are currently possible. 
Tf the searcher makes a CZ detection, KCZO = -l. 
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ROMR/inba@ger/..... The number of the most distant CZ 
in which the target can currently make a detection. 
If KCZT = 0, no CZ detections are currently possible. 
Lr tHe target makes a CZ detection, KOZT = -1. 


RR/*/int@o@r/..... The number of the searcher track 
anchor point at which the search pattern begins to 
repeat. This is the "return point" from which 
subsequent patterns will start after the first full 
pattern 1S executed (input). RESTR: positive, less 
than NP. See NEXT(I) and NP. 


mo/* /integer/.=... The number of range/propagation 
loss data points used to define the searcher direct- 
path (DP) propagation loss curve (input). 

RESTR: greater than 1, and less than 21. 


MODE/integer/..... Iné@eator variable. If MODE = l, 
Searcher 1S at drift speed. If MODE = 2, searcher 
1S at Sprint speed. 


MSELO/integer/..... The MSEL counter for the searcher. 
MSELT/integer/..... The MSEL counter for the target. 
Mmp/*/integéer/..... The number of range/propagation 


loss data points used to define the target direct 
path-propagation loss curve (input). 
RESTR: greater than 1, and less than 2l. 


MBWMNS /*/integer/..... Thenumber of sections, each 
with a fixed probability, that define the initial 
target distribution across the top of a barrier choke 


point (input). RESTR: greater than 1, and less 

tiene 21 . 

NBOTH/integer/..... The number of simultaneous 
detections. 

NCZDO(I) /integer/..... In the main program, the number 
of searcher detections made in the Ith convergence 
zone. Converted in subroutine SINKEM to the fraction 


of searcher detections in the Ith CZ. 


NCZpe (i) /integer/..... In the main program, the number 
of target counter-detections made in the Ith conver- 
gence zone. Converted in subroutine SINKEM to the 
fraction of target detections in the Ith CZ. 
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MCAO/*/ damtc@iger/...%.% The number of searcher convergence 


zones (input). RESTR: non-negative, less than 6. 
WOZT/?T/ismftegéer/..... The number of target convergence 
zones (input). RESTR: non-negative, less than 6. 
NDO/integer/..... The number of detections made by 
searcher. 

NDT/integer/..... The number of detections made by the 
Garget. 

Weer (lL) /integer/..... The next anchor point in the 
searcher track after point I. The program sets 

Meer (lt) = I+] for I= 1 to NP-1, and NEXT(NP) = KP. 


Thus, after one complete execution of the search 
pattern, the searcher transits from point NP to KP, 
and repeats the pattern starting at KP. 


NL/integer/..... Pointer to the next searcher track 
leg (similar to NLEG). 


MBEG/imbeger/..... A pointer to indicate the searcher 
position in the search pattern. The searcher is 
always between anchor point NLEG and point NEXT(NLEG). 


me /*/mnteger/..... The number of relative bearing/ 
figure-of-merit data points used to define the radial 
distribution of target Source level resulting ina 
radial distribution of searcher figure-of-merit 


(Input). RESTR: greater than 2, less than 51. 
NONE/integer/..... The number of replications in which 
neither ship makes a detection prior to TMAX. 

fee, *7integer/..... The number of searcher track anchor 
Beints (input). RESTR: greater than 1, less than 
Se. 

NeePp7*/interger/..... The number of replications to be run 
fnput). RESTR: positive, less than 5001. 

NWTDO(L) /real/hours/..... The time of the Ith Searcher 
detection. 

Pewr(l) /real/nours/..... The time of the Ith target 
counter-detection. 

OFOM(I)/real/decibels/..... Through the EQUIVALENCE 
Statement, OFOM(1) = FOMOD, OFOM(2) = FOMOS. The 


subscript of this variable in the program is MODE. 
See MODE, FOMOD, FOMOS. 
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Olb( 3.) /* / read /decibeds/...«: Direct-path propagation 
loss at range RO(I) for the Searcher (input) . 
RESTR: OL(I) must increase as range increases. 
Git) < OL(I#]) . 


PREN (1) /*/medlL/..... The probability of target inclu- 
Buen an the Jth geection when the user defines inztial 
ia@eral target distrabutition, I = l,...,NBINS, (imput) . 
RESTR: non-negative, less than or equal to 1.0. 

Sum of all PBEN(I) = 1.0. 

Rezo (LJ) /*/ Beal /navtwcal miles/ ... «« Range to the 
inner (J=l1) and outer (J=2) edges of the Ith conver- 
gence zone for the searcher, I=l,...,NCZO, (input). 


Por R: RCZO(H,1) < RCZO(I,2) < RCZO(I+1,1) < 
REZO(I+1,2),. etc. 


RCZT(1I,J) /*/real/nautical miles/..... Range to the 
inner (J=l1) and outer (J=2) edges of the Ith conver- 
gence zone forthe target, I = 1,...,NCZT, (input). 


Peewee em () 6 << REAT(ET,2) < RCZT(I+1,1) < 
Meee it+l,2), etc. 


eee t/loOgucal/..... RFLCT is FALSE unless a target 
reflection takes place. 

RMAX/real/nautical miles/..... The maximum of RNGO,RNGT. 
RNG/real/nautical miles/..... Current range from 


Searcher to target. 


RNGO/real/nautical miles/..... Searcher current direct 
path detection range. 


RNGT/real/nautical miles/..... Target current direct 
path detection range. 


RNTDO(I) /real/nautical miles/..... Range between 
Searcher and target when searcher detects target. 


RNTDT (I) /real/nautical miles/..... Range between 
Searcher and target when target detects searcher. 


RO(1I)/*/real/nautical miles/..... Range associated 
with searcher propagation loss OL(I) (input). 
RESTR: Must be entered in ascending order, RO(I)< 
RO(I+1). 
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Re (1) /7*/veal/nautical mirLea/..... Range associated 
with target propagation loss TL(I) (input). 

RESTR: Must be entered in ascending order, RT(I) < 
RT(I+a). 


RTCC/*/real/inverse hours/..... Rate of target course 
changes. 1/RTCC is the mean of an exponential distri- 
bution from which the time interval to target course 
change is drawn (input). RESTR: positive. 


RTSC/*/real/inverse hours/..... Rate of target speed 
changes. 1/RTSC is the mean of an exponential distri- 
bution from which the time interval to target speed 
change is drawn (input). RESTR: positive. 


SRED/*/integer/..... Psuedo-random number generator 
seed (input). RESTR: non-negative. 


SIGMA(I)/*/real/decibels/..... Scale parameter for the 
th “acoustic Fluctwation process, I=1,2,3. 
RESTR: non-negative. 


Sankx/realy/..... The y-component of a unit vector in 
the direction of target motion. 


BO Treat/KnoOts/y..... Current searcher speed. 


S@p/*/real/knots/..... Searcher drift speed (input). 
For constant speed search, set drift speed equal to 
Sprint speed (SOD = SOS). 


SOS /*/real/knots/..... Searcher sprint speed (input). 
For constant speed search, set sprint speed equal to 
drift speed (SOS = SOD). 


were logicaly.....- This variable is TRUE if the user 
has specified a lateral target distribution other 
than uniform across the barrier (barrier scenario 


only). Otherwise, it 1s FALSE. 
SPEED(I) /real/knots/..... Through the EQUIVALENCE 
statement, SPEED(1) = SOD, SPEED(2) = SOS. The 


subscript of this variable in the program is MODE. 
See MODE, SOS, SOD. 


ST/real/knots/..... Target current speed, drawn from 
a uniform distribution on the interval (STMIN,STMAX). 
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Slaeenl , STACK 246, STACKS ,STA@KG/ readl/ sia. Beach 39 @ stack 
of 5000 psuedo-random numbers which are sequentially 
dereferenced for program use by the stack counters 
CiyG2,¢3,Cé. SPACRL ve distributed exponentially, 
mean = 1.0. STACK2 is sistributed normally, mean = 
0.0, variance = 1.0. STACK3 is distributed unifornm- 
aliy on the imterval (0,1). STACKd is distributed 

as specified by the user. 


START/*/ real /hours/....: The latest time-late on 
barrier time. The time-late on barrier for the 
searcher is distributed uniformly on the interval 
(0,START), with an assumption that the target travels 
at mean speed from time zero until the searcher is on 
the barrier (input). RESTR: non-negative. 


Sueemo L/logmuecakb/..... This variable is FALSE if no 
statistical analysis is to be done. Otherwise it 
Ve TREE . 


Cie 2/ bhogread/ ..a.. This variable is TRUE if section- 
ing of the data 1S to be accomplished. Otherwise, 
it is FALSE. 


STINC/real/knots/..... The interval of speed which 
defines the uniform distribution of allowable target 
speeds (STMAX - STMIN). 


STMAX/*/real/knots/..... Target maximum speed (input). 
RESTR: STMAX must be greater than or equal to STMIN. 


STMIN/*/real/knots/.....Target minimum speed (input). 
RESTR: STMIN must be less than or equal to STMAX. 


TBIG/real/nautical miles or hours/..... TBIG is used 
to compute the simulation time step. Initially, TBIG 
is a range = RNG - RMAX, and is positive if, and only 
if, no direct path detection is possible. If no 
detection is possible (CZ or DP), TBIG is set to the 
minimum closing distance such that either a detection 
or counter-detection 1s just possible (DP or C2) 
based on the current acoustic conditions. TBIG is 
converted to a time interval by dividing by the 
current combined speeds of the searcher and target, 
which assumes a worst case head-to-head closing 
Situation. The simulation time-step is then set to 
the maximum of TBIG or 0.05 hour. 


NEG teal /nourey..... Time of next searcher course 
change. 
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TD) /)”/ me @l/hemmgs / «. ..« Searcher drift time (input). 
FESTR> positive. 


TFOM(I) /real/decibels/..... Through the EQUIVALENCE 
Statement, TFOM(1) = FOMTD, TFOM(2) = FOMTS. The 

Subscript of this variable in the program is MODE. 
see MODE, FOMTD, FOMTS. 


Tae T/ real /madwans/..... Target psuedo-course. THETA 
is measured from the x-axis in a counter-clockwise 
di neceon. 


TIFPL(1) /peal/hours/.«..%* Ber Meal ,2,3, TIFL(L) ws the 
time at which AFL(I) changes in the LSJ process, or 
the next sample time for the GMA process. For 

I = 4-7, through the EQUIVALENCE statement, 

Prec) * TCC, TIFPL(S) = TTCC, TIFL(6) = TTSC, 
el) = TSC. 

Tne (lL) /real/hours/..... Through the EQUIVALENCE 
statement, TIME(1) = TD, TIME(2) = TS. The subscript 


of this variable in the program is MODE. See MODE, 


TINC/real/hours/..... The time at which the next 
Simulation time-step ends. 


mmect)/*=/real/decibels/..... Propagation loss at range 
RT(I) for the target (input). MTL(I) must increase 
with increasing range. TL(I) must be less than 

me (t +p ae). 


THAST/real/hours/.... The time at which the 
Gauss~Markov error process was last evaluated. 


TMAX/*/real/hours /. . ss» Maximum allowed search time 
per replication. 


TNOW/real/hours/..... Current time in a replication. 
wel) /real/nours/....: The time at which the figure- 
of-merit sample for the searcher, FO(I), was taken. 
TOXMAX/real/nautical miles/..... Twice the value of 

XMAX. 

TOYMAX/real/nautical miles/..... Twice the value of 

YMAX. 

TS/*/real/hours/.....Searcher sprint time (input). 


REotR: positive. 
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tTeCc/real /n@mmes/:i...s% Time of next searcher speed 
change. 


or l)/ee@al/noursy..../fhe time @t whileéh the tafget 
figure-of-merit sample, FT(I), was taken. 
TFCC/ue@al./hours/.:i... Time of next target course 
change. 

TTSC/real/hours/..... Time of next target speed 
change. 

TwOePT/real/..... Twice the value of T. 

UR /meal/knots/..... Tre x=component of current taxget 


velocity vector. 


Uy real /kMOte/..... The y-component of current target 
velocity vector. 


VX/real/knots/..... The x-component of current 
Searcher velocity vector. 


VY¥/real/knots/..... The y-component of current 
searcher velocity vector. 


RBEN() /*/real/nautical miles/..... The distance from 
the origin (x = 0.0) to the right-most x-value of 
the section containing target probability mass 
PBIN(I). RESTR: positive, XBIN(NBINS) = XMAX. 


XMAX/*/real/nautical miles/..... Length of the search 
area rectangle in the x-direction (input). 
RESTR: positive. 


XO/real/nautical miles/..... The current searcher 
x~position. 


XODT(I)/real/nautical miles/..... The x-position of 
the searcher when the searcher makes the Ith 
detection of the target. 


XOT/real/nautical miles/..... The projected searcher 
position based on the current position (XO) and the 
computed simulation time-step. 


XOTD(I) /real/nautical miles/..... The x-position of 


the searcher when the target makes the Ith counter- 
detection of the searcher. 
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MP (1) /*/real/namtical miles/..... The x-position of 
the Ith searcher track anchor point (input). 
RESTR: non-negative. 


MT / real /newtueced miheg/.«... The current x-position 
of the target. 


XTDO (CI) /real/namtacal, mi 16S / . nom. The x-position of 
the target when the target makes the Ith counter- 
detection of the searcher. 


PIOD(I) /real/navwiical miges/..... The x-poSition of 
the target when the searcher makes the Ith detection 
of the target. 


YMAX/*/real/nautical miles/..... Length of the search 
area in the y-direction (input). 

RESTR: positive. 

YO/real/nautical miles/..... Current searcher 


v= pesn. tien. 


YODT(I) /real/nautical miles/..... The y-position of 
the searcher when the searcher makes the Ith 
detection of the target. 


YOT/real/nautical miles/..... Projected searcher 
y-position based on the current position (YO) and 
the simulation time-step. 


YOTD(I) /real/nautical miles/..... The y-position of 
the searcher when the target makes the Ith counter- 
detection of the searcher. 


YP(I)/*/real/nautical miies/..... The y-position of 
the Ith searcher track anchor point (input). 
RESTR: positive. 

YT/real/nautical miles/..... Current y-position of 


the target. 


YEwe(l) /reas /nautical miles/«.... The y-position of 
the target when the target makes the Ith counter- 
detection of the Searcher. 


YTDO(I) /real/nautical miles/..... The y-position of 


the target when the searcher makes the Ith detection 
of the target. 
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ABBREVIATIONS AND ACRONYMS 


This appendix contains an alphebetical listing of abbrev- 
lations and acronyms uSed in this thesis. 
1. ASW: Antisubmarine Warfare. 
meeecewv(X,Y): If X,Y awe randem variables, this is the 
covariance function. If X,Y are realizations of an 


autoregresSive time series, this is the autocovari- 
ance function. 





3. CPU: Central Processing Unit. 


4. CZ: Convergence Zone. Usually used to describe the 
propagation mode of an acoustic Signal. 


S. db: Decibel. 


| 6. DP: Direct Path. Usually used to describe the propa- 
| gation mode of an acoustic signal. 


7. §E(X): The expected value of the random variable X. 
ee PR: Exchange Ratio. 


9. EXP(x): An exponential distribution with rate 
parameter x. 


10. FOM: Figure-of-Merit. 


ll. IMSL: International Mathematics and Statistics 


Library. 
12. Le: Environmental noise. 
13. Ls: Target noise source level. 


14. MOEs: Measure Of Effectiveness for the Searcher 
15. MOEt: Measure Of Effectiveness for the Target. 


16. MSEL: Minimum Signal Excess Logic. 
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N(x,y): A normal distribution with mean x and 
variance y. 


NONIMSL: Non-International Mathematics and Statistics 
Library. 


NPGS (NPS): Naval Postgraduate School. 
Nrd: Recognition differential. 
PASS: Passive Acoustic Search Simulation. 


PCD: Probability of Counter-Detection. 


PCDCZ: Probability of Counter-Detection in a 
Convergence Zone. 


PCDCZI: Probability of Counter-Detection in the Ith 
cConm@rgence zone. 


PCDDP: Probability of Counter-Detection by Direct 
Path propagation. 


PD: Probability of Detection. 
PDCZ: Probability of Detection in a Convergence Zone. 


PDCZI: Probability of Detection in the Ith. 
Convergence Zone. 


PDDP: Probability of Detection by Direct Path 
propagation. 


0(X,Y): If X,Y are random variables, op is the 
linear correlation coefficient. If X,Y are 
realizations of an autoregresSive time series, 

then 9 is the autocorrelation function. 

Rs/Rt: The mean detection/counter-detection range. 


(R50)s/(R50)t: The median detection/counter- 
detection range. 


SE: Signal Excess. 
Ts/Tt: Mean time to detection/counter-detection. 


(T50)s/(T50)t: Median time to detection/counter- 
detection. 
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TMA: Target Motion Analysis. 


Wi(@e.pb:). 


Var(X): 


A uniform distribution on the interval 


The variance of the random variable xX. 
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APPENDIX D 


APPLICATION OF THE PASSIVE SONAR EQUATION TO THE PASS MODEL 


The pasSive sonar equation may be written as: 
SNR = LS - Nw - Le 
where: 


1. SNR = the signal-to-noice ration at the processor 
Sutpw@t (db). 


2. Ls = the target source level (db). 


3. Nw the signal propagation loss (db). 


4. Le = the background noise compensated for by sonar 
aqerectavity (db). 


Now, make the following definitions: 

1. Nrd = RECOGNITION DIFFERENTIAL (db). Nrd is that SNR 
required for a probability of detection of 0.5. It 
can be thought of as the SNR required for a target to 
be recognized as a valid contact 50% of the time. Nrd 
1s generally a function of sonar processor performance. 

2. FOM = FIGURE-OF-MERIT (db). FOM is that amount of 
Signal propagation loss (Nw) that results in a 
probability of detection of 0.5. 


3. SE = SIGNAL EXCESS (db). SE is the algebraic differ- 
ence of FOM and Nw (SE = FOM - Nw). 


If we now substitute Nrd for SNR and FOM for Nw in the above 


equation, and solve for FOM, we have: 


FOM = Ls - Le - Nrd 


Thus, the result 1S: 


SE = 0 <= Nw = FOM=> Pd 
SE > 0 = Nw < FOM™ Pd > 
SEe< 04> Hw > FOM<= Pd < 0 
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In PASS, we obtain a certain detection whenever the sum 
of the deterministic Signal excess and the error term is 
greater than or.equal to zero. The error term has a zero 
mean, and the probability that 1t is greater than zero is 
0.5. Thus, when the deterministic Signal excess iS zero, the 
probability of detection in PASS is 0.5, which 1s consistent 


with the passive sonar equation results above. 
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APPENDIX E 


CerewreR SPHULATION OF THE STOCHASTIC ERROR PROCESSES 


In this appendix we examine the modeling of the acoustic 
Fluctuation process independent of all other events in the 
Simulation. In effect, we will ignore all events which 
affect the simulation time-step except for acoustic 


fluctuations. 


ae 0U ee, LAMBDA-SIGMA JUMP PROCESS 
In PASS there are three fluctuation processes simultane- 


ously in progess. 


X(t) = searcher local error process. 
Y(t) = global error process 
Z(t) = target local error process 


Each process is defined by a rate parameter (Ay 5 hos h3) 
and a scale parameter (85, O45, 03). For illustrative 
purposes, assume we initially start the process at time = T 
with independent draws from normal distributions parameter- 


ized by the scale parameters. That is: 


- 2 

X(T) = Dies where Gi, uN (0,0,) 
2 

imc) = where s N (0,07) 
cc) Ny ny ( 
mot) = ie where n, VN (0,03) 


Bs)is 








Also, at time = T, we determine the time interval to the 
next fluctuation level change, s, by independent draws from 
exponential distributions parameterized by the rate 


parameters. That is: 


Sy a where i. EXP (4) 


WN 
II 


s leg where Ty \ EXP (i.) 


5. = i where tere EXP () 3) 


Me actual time of the fluctuation changes would then be: 


fear, * = 
aoe 5S ot T 
y y 

es 2a * fT 


If no detection takes place at time = T, the next time a 
detection can take place is the time of the first fluctuation 


change. Therefore, the simulation time=step , At, would be: 


yiSz) 


At = min {S.,S 
All simulation events (e.g. platform motion) would then be 
executed with time-step At, and the simulation time would be 
incremented by At. 


iMye Llwectuation level associated with At is then 


changed. Suppose, for example, that t.<ty st: Then, we 
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would change the X fluctuation level by replacing the 
previous level with a draw from N(O,05), and obtain a new 
time interval to the next change of X by making a draw from 
EXP(A,), from which the time of the next change in X 1s 
calculated. If no detection takes place, a new time-step is 


computed as before, and the process 1S repeated. 


B. THE GAUSS-MARKOV PROCESS 
Unlike the LSJ process, the GMA process involves continu- 
ous sample paths. The defining function of the GMA process 


2S Of the form: 


an X(t) + g(s)n, 


X(t + s) 


where: 
i. Mise is a normally distributed random variable with 
zero mean and variance of 
2. g(s) is a function of the time increment, s, which 
has the properties: 


lim g(s) = 0 and lim g(s) = l 
S+O S +00 





we know that: 
Var[X(t)] =o 
SO 


AS 


Var[x(tt+s)] = of = Var [fe Rey + g(s)n,] 


6 








= xS 
Since e and g(s) are not random variables and X(t) and a. 


are independent, we can write: 


of = aes Ver iese) ] + [g(s)]° var[n,] 
a4 = anes of + [g(s)]* of 


ee (j-e7 248) 1/2 


muse, the functional form of GMA is: 


X(t+s) = oX(t) + aie a is 


where: 


The difficulty in using this process in a computer simu- 
lation is deciding on the time-step. The more accurately one 
wants to model the process, the smaller the time-step 
required, but small time-steps cause excessively long run 
times. As a practical compromise, the PASS model uses the 
identical method for determining fluctuation time-steps in 
both the LSJ and GMA process (i.e. the GMA process is 
sampled at exponentially distributed random times governed 
by the rate parameters). The difference in applying this to 
the GMA process is that each time one fluctuation process is 


changed in LSJ, all three are updated in GMA. 


vo 








For example, in the previous section where X was changed, 
the GMA process would update X, Y, and Z based on a time 


, 1 


increment Sur and independent draws to Nye iv 


me 
C. GENERATION OF PSUEDO-RANDOM NUMBERS 

Psuedo random numbers used throughout PASS are generated 
using the NPGS random number package LLRANDOMII. This 
package was used, vice a portable package (e.g. International 
Mathematics and Statistics Library) because of the improve- 
ment in speed in number generation. The IBM Assembly 
Language LERANDOMII routines are between three and eight 


times as fast as the FORTRAN IMSL subroutines. 


Pee G@eRSGATION OF THE ERROR PROCESSES 
Suppose we define two compound error processes, A and B, 


as follows: 


PG ) X(t) + Y(t) 


B(t) fate) + Z(t) 


where X, Y, Z are LSJ or GMA processes. Dropping the "t" 
notation for convenience, a computational formula for the 


covariance is: 


Cov(A,B) = E[A°*B]-E[A] -E[B] 


II 
© 


mmcestmmce E[A] = E[B] 


EB [XY+Y°4XZ+YZ] 


Cov(A,B) = E[A°-B] 
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GevuwB) = Cov(xX,y¥) + a + Cov. (pe) & Gov (FZ) 


Geodesy) = a7 
y 
Since X,Y,Z are independent, and E[yY] = 0. 


The linear coefficient of correlation is defined as: 


eS) = Cov (A,B) 
AB 
where 
= 2 Deeliy/ 2 
- 2 relay 2 
- = (a5 - d3) 


Peyeeme coefficient of correlation for the compound error 





process is: 


= 1 5 
2 ’2 
Woete that: 1. lim o(A,B) = 1 
O O 
=~ = 2 
2 2 


which satisfies the intuitive and practical requirement that 
as the global fluctuation becomes dominant, the error signal 


experienced by both sensors will be highly correlated. 


De. lim o(A,B) = 0 
S aeaeene 
— , —— =O 
27 “9 


which satisfies the intuitive and practical requirement that 
as the local fluctuations dominate, the error signal experi- 


enced by the sensors will be almost uncorrelated. 
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APPENDI. ? 


A WP REB=OUT=O0OF-FIVE DETSECTION CRETERTA. ODEL 


The Minimum Signal Excess Logic (MSEL) starts in state 
zero (MSEL 0). Each time the acoustic environment is 
samples, the MSEL state 1S incremented by one if detection 
marmeessible (SE > 0). If detection is not poSsSible (SE < 0), 
then the MSEL counter is decremented by one, unless MSEL is 
in State zero. The lowest MSEL State iS zero, and detection 
occurs the first time the MSEL state reaches three. This 
logic reSultS in a detection occuring if, and only if, three 
of the last five samples of the acouStic environment had 
cer2 0. 

MSEL is equivalent to the Markov chain shown in Figure 
F.1, where the transition probabilities are, in general, a 
function of time (i.e. they will change each time the 
acoustic signal is sampled). In uSing the LSJ process, this 
presents no problem in that the fluctuation process 1s 
Static between samples. In uSing the GMA process, however, 
the fluctuation process iS dynamic between sample times, and 
therefore the rate of Sampling will have an effect on the 


time of detection. The PASS model accounts for this in two 


ways: 
1. The rate of environment sampling is the same for LSJ 
and GMA. See Appendix E. 
2. Whenever either searcher or target are not in MSEL 0, 
the maximum PASS time-step defaults to 0.05 hour (3 
minutes). 
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For analysis purposes, assume all the transition prob- 


@eelities are constant and equal to p. bet q = (l-p) and 
tet Si = the mean number of time increments (samples of the 
acoustic environment) to go from state i to state 3. Reach- 


ing the absorbing state, state 3, is equivalent to detection. 


We can write a series of balance equations as: 


”n 
II 


5 p(1l) + q(l + S,) =] + qs, 


— 
II 


1 p(l + pum o(1 + So) =e + pS, + qSo 


Se p (1 S,) ae ( loet Se) 1 + pS, + qSo 


If we solve this system of equations for So, we have: 


So = 3(2 + =5) 
P 


Since So is the mean number of time steps to go from 
MSEL state zero to MSEL state 3 (for a fixed p), we can 
place an upper bound on the expected time to go from MSEL 0 


to MSEL 3 as: 
To = 3So (in minutes) 


Since the maximum time-step is PASS is three minutes when 
the searcher or target MSEL counter 1s greater than zero. 
Table VIII shows representative values of So and To as a 
function of the instantaneous probability of detection. 
Keeping in mind the simplifying assumptions made, the values 


of To are analogous to maximum expected integration times. 
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The MSEL Model as a Markov Chain. 


TABLE VIII 


MSEL MAXIMUM EXPECTED INTEGRATION TIMES 


J Mm & WwW WW 


he 
208 
43. 
iS > 
O20. 


50 


. 000 
94 
-453 
oi tS 
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000 
625 
704 
000 
000 
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APPENDIX G 


CODE FOR EXECUTIVE PROGRAM TO RUN PASS 


& TRACE 

CP TERMINAL LINESIZE 80 

ERASE PASS OUTPUT 

ERASE PASS HOLD 

ERASE PASPLT DATA 

COPYFILE PASS DATA Al = HOLD = 

ERASE PASS DATA 

FILEDEF 06 DISK PASS OUTPUT 

FILEDEF 07 DISK PASS DATA 

FILEDEF 08 TERM 

FILEDEF 09 DISK PASS HOLD 

FILEDEF 10 DISK PASPLT DATA 

CLRSCRN 

&TYPE PASSIVE ACOUSTIC SEARCH SIMULATION *** “Pass" 

& TYPE 

&TYPE EXECUTING PROGRAM "PASS" ** "PASS" MUST BE FORTRAN-G 
OR FORTRAN-H 

&TYPE PRECOMPILED ON YOUR A-DISK. THE NONIMSL LIBRARY MUST 
BE AVAILABLE. 

&TYPE 

&TYPE DATA FILE MANAGEMENT: 

& TYPE . 

&TYPE A. INPUT DATA FILE "PASS DATA": 

&TYPE THIS FILE CONTAINS THE MOST RECENT DATA FOR RUNNING PASS. 

&TYPE IT IS CREATED, OR MODIFIED, AFTER RUNNING THE INTERACTIVE 

&TYPE DATA INPUT SECTION OF PASS. SEE FURTHER NOTES ON OUTPUT 
DATA 

&TYPE FILE "PASS HOLD". 

& TYPE 

&TYPE B. OUTPUT DATA FILE "PASS HOLD": 

&TYPE UPON EXECUTIO OF PASS, IF FILE "PASS DATA" ALREADY 
EXISTS, 

&TYPE IT IS COPIED INTO “PASS HOLD" UNCHANGED. THUS IF YOU 
CHANGE THE 

&TYPE INPUT DATA DURING THE INTERACTIVE SECTION, THE CHANGED 
DATA WILL 

&TYPE BE IN "PASS DATA", AND THE UNCHANGED DATA WILL BE IN 
"PASS HOLD". 

&TYPE 

&TYPE C. OUTPUT DATA FILE “PASPLT DATA": 

&TYPE THIS FILE CONTAINS ALL THE DATA, IN THE PROPER FORMAT, 
NEEDED 

&TYPE TO GENERATE GRAPHICAL OUTPUT USING PROGRAM "PASSPLT", 
AND THE 


oS 








&TYPE DISSPLA GRAPHICS SYSTEM. 
&TYPE 
fexre D. OUTPUT DATA FILE “PASS OUTPUT” : 
Stare THis FILE CONTAINS THE STATISTICAL OUTPUT FROM THE 
MOST RECENT 
STYPE RUN OF PASS. 
&TYPE 
a a a i Se 
LOAD PASS (START 
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APPENDIX H 


FORTRAN CODE FOR PASS 
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